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We study the equation of motion for the Noether current in an electron gas within the 
framework of the Schwinger-Keldysh Closed-Time-Path formalism. The equation is shown to 
be highly non-linear and irreversible even for a non-interacting, ideal gas of electrons at non¬ 
zero density. We truncate the linearised equation of motion, written as the Laurent series in 
Fourier space, so that the resulting expressions are local in time, both at zero and at small 
finite temperatures. Furthermore, we show that the one-loop Coulomb interactions only 
alter the physical picture quantitatively, while preserving the characteristics of the dynamics 
that the electric current exhibits in the absence of interactions. As a result of the composite 
nature of the Noether current, composite sound waves are found to be the dominant IR 
collective excitations at length scales between the inverse Fermi momentum and the mean 
free path that would exist in an interacting electron gas. We also discuss the difference and 
the transition between the hydrodynamical regime of an ideal gas, defined in this work, and 
the hydrodynamical regime in phenomenological hydrodynamics, which is normally used for 
the description of interacting gases. 
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I. INTRODUCTION AND MOTIVATION 

Hydrodynamical description of many-body systems with collective excitations is a well-established 
and widely applied method, which successfully combines microscopic and phenomenological in¬ 
formation about the physical system. Despite the numerous successes of this mixed scheme, 
the systematic, microscopic derivation of the hydrodynamical equations is not fully understood. 
Rather than deriving them directly, one usually calculates transport coefficients from the micro¬ 
scopic description and uses their values in the equations of phenomenological hydrodynamics, i.e. 
the gradient expanded constitutive relations [1-3] . 

Much recent research has been focused on understanding hydrodynamics in terms of an effective 
theory of Goldstone modes [4-8], including works that proposed ways of incorporating dissipation 
into the effective description [9-12]. However, none of these approaches provides us with a detailed 
microscopic view of the origin of the effective field theory of hydrodynamics. This present work grew 
from a desire to bridge this gap and to gain new insight into the ab initio elements of a hydrodynamic 
system, derived directly from a microscopic quantum field theory, including dissipation. In essence, 
this paper is one of the simplest examples of a microscopically derived effective Schwinger-Keldysh 
field theory with the structure of the classical theory explored in [9]. 

An important simplifying feature of this work is that our present goal is to understand the role of 
many-body correlations that characterise the hydrodynamic flow, independently of the interactions 
among the participating particles. The suspicion that highly non-trivial correlations might be 
present even in the absence of interactions comes from the simple fact that gauge-invariant (i.e. 
charge conserving) observables are composite operators, and hence the structure of their connected 
Green’s functions is richer than for those of elementary operators. We will argue that similar types 
of phenomena are not restricted only to quantum systems but that they are present also in classical 
mechanics. 

Our main goal in this paper is to find the equation of motion for the expectation value of the 
current operator, (j^), in an ideal gas of electrons at non-zero density, both at vanishing and low 
temperatures. In the absence of any phenomenological input, the relevant effective theory must 
be such that (j^) satisfies its variational equation. It is the calculation of expectation values that 
requires us to go beyond the conventional formalism of quantum field theory and perform the 
calculation in the Schwinger-Keldysh Closed-Time-Path (CTP) formalism [13-17]. 

The CTP effective action will be calculated within the spirit of the Landau-Ginzburg double 
expansion, by organising each order of the expansion in powers of the amplitude. This expansion 
gives rise to a Laurent series and orders the terms in the equation of motion in increasing powers of 
the wave vector. This should be compared with the phenomenological approach to classical hydro¬ 
dynamics where the equations, expressing the conservation of the energy-momentum and charges, 
are also constructed by using the traditional double expansion [1]. These two approaches are fully 
analogous with the exception of three important differences. Firstly, the equations derived from 
microscopic physics are closed without the requirement of additional thermodynamical considera¬ 
tions. Local equilibrium in an infinite, homogeneous system is automatically ensured so long as the 
perturbations are weak and slowly varying. As a result, the usual difficulties of the phenomenolog- 
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ical approach to zero temperature hydrodynamics appear as a problem in performing the Legendre 
transform in deriving the effective action. Secondly, a new element of this approach is the fact that 
the derivation reveals a spatial non-locality induced by the Fermi surface, which is reflected by the 
appearance of the negative powers of the wave vector, l/|q|, in the Laurent series. Such terms, 
generated by loop-integrals, are absent in the usual phenomenological approaches, which are based 
on the naive Taylor series expansion of classical physics, i.e. in oj and q} The third difference 
concerns the definition of the hydrodynamical regime. In the phenomenological approach, the 
relevant regime is defined at length scales where the local density, current and thermodynamical 
potentials can be defined. In our approach, based solely on local expectation values, there is a well 
defined effective action, valid at any length scale, well beyond the minimal length, set by the UV 
cutoff of the underlying quantum field theory model. We define the hydrodynamical limit by the 
requirement that the equations of motion be local in time and display integro-differential structure 
in space. 

The transport coefficients, entering into the phenomenological hydrodynamical equations, can 
be derived from quantum field theory by using Kubo’s formulae [3, 21, 22]. Our generating func¬ 
tional for the connected CTP Green’s functions of the current, calculated at the quadratic order, 
reproduces the usual Kubo formulae. It is then the next step in our work that goes beyond the tra¬ 
ditional treatment of the electron gas. We perform the functional Legendre transform and calculate 
the effective action for the current, which is what produces closed equations of motion without the 
need for any thermodynamical input. Even though the derivation of the standard Kubo formula 
for electrical resistivity [21] is the first step in obtaining the dynamical equation for (j^), this line 
of research has not been pursued in a systematic manner to our knowledge. 

The behaviour of a Fermi liquid is normally formulated by Landau’s phenomenological the¬ 
ory [23-25], which is based on the idea of quasi-particles. However, the charged, local fields that 
represent the quasi-particle excitations around the Fermi surface are not Hermitian physical observ¬ 
ables. This problem is circumvented in Landau’s argument by carefully constructing the effective 
dynamics near local equilibrium in terms of the scattering processes of the quasi-particles. The 
complication in deriving such a theory from microscopic physics in a systematic manner is the 
necessity for using non-local transition amplitudes, which are given by the gauge invariant residues 
of the amputated, connected Green’s functions. Instead, we wish to formulate an effective theory 
in terms of local observables to understand how phenomenological equations of hydrodynamics 
arise. We must therefore rely on neutral, composite operators. The obvious choice, the current 
j^{x) = makes our effective theory formally different from Landau’s Fermi liquid 

theory. 

The main observation we wish to present in this work is that a non-interacting, ideal gas of 
electrons displays complicated correlations and dynamics when one analyses the behaviour of the 
electric current. It exhibits features that are usually characteristic of interacting systems, such as 
the presence of a sound mode. To better understand the role of interactions, imagine a calculation 
of the AG electrical conductivity in an electron gas. This can be done by introducing an external 

^ An example of the complexity of hydrodynamics beyond an analytic Taylor expansion in uj and q, which has been 
well known and also arises in the computation of current-current correlation functions (away from large-A^) are 
the so-called “long-time tails”, which exhibit non-analytic behaviour in to [3, 18-20] . 
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plane-wave electric field with the amplitude E = {E, 0, 0) and then calculating the response of the 
system, namely the amplitude of the induced current j = (j, 0,0). Perturbatively, the result has 
the form 


j 


OO OO 

— C ^ ^ ^ ^ Cmn 
m=l n=0 


eE 



( 1 ) 


expanded in terms of two small parameters, the QED coupling strength, e, and the ratio of the 
external perturbation, eE, to the Fermi energy, ep. The double expansion reflects the double 
role the electromagnetic interactions play in this problem; on the one hand, e is the fundamental 
interaction among the charges in the gas and on the other hand, it is the interaction used to 
diagnose the gas. The key point is that the dynamics of an out-of-equilibrium gas behaves differently 
depending on the order the two independent e —)■ 0 and eE/ep —^ 0 limits are taken. 

To characterise the gas itself, without the large interference of the external probe, one goes into 
the limit of the linear-response regime. Namely, one carries out the limit eE/ep —)■ 0 first in such 
a way that the transport coefficient, i.e. the conductivity, becomes well defined and finite. 





OO 


n=0 



( 2 ) 


In this limit, the interactions between the gas and the observational probe are severely simpli¬ 
fied. However, the interactions among the charges still play an important and complicated role in 
determining the coefficients cin- To simplify the contribution from the fundamental interactions, 
one usually considers the un-physical limit in which the QED interaction becomes extremely weak. 
The leading-order result, aje^ ~ cio/ep, is then sufficient for describing the dynamics of a (nearly) 
ideal gas in the presence of a weak external perturbation eE. 

Alternatively, one may be interested in the way the elementary interactions of QED dress the 
non-interacting, ideal gas in the presence of external perturbations. This scenario corresponds to 
the other order of limits where e —)• 0 is performed for finite eE. In this approach, the external 
perturbation has to be restricted to an AC form in order to avoid the instabilities of non-interacting 
particles subject to homogeneous external forces. This is the strategy followed in this paper in which 
we derive the linearised equation of motion for the electric current in the presence of an external 
perturbation, within the one-loop approximation in QED. The external source, a^(x), is fictitious 
in our scheme. It is used to drive the electron gas from its equilibrium state at U = —oo to a 
desired, non-trivial state at t = 0 when the source is switched off, i.e. a^(x) = 0 for t > 0. The 
equation of motion, satisfied in the absence of the source, is therefore well defined so long as it is 
local in time. Since it is only possible to find a local equation of motion for sufficiently slow motion 
of the gas, we will restrict our attention to external sources with slow time dependence. 

The non-trivial correlations in a non-interacting, ideal gas appear due to the redistribution of 
the energy-momentum into its non-interacting normal modes, when the dynamics is diagnosed 
by composite operators. The initial energy-momentum is injected into the gas from the external 
source. Such a redistribution generates new collective modes with non-trivial dispersion relations, 
leading to new composite sound waves. Dissipation manifests itself in the spread of the flow pattern 
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in space, as a function of time. As a result, this type of a sound wave is different from both the 
zero and first sound modes of the usual Fermi liquid. 

There are three major parts of this paper: the discussion of inertial forces that arise from non¬ 
linear coordinate transformations, the argument that effective theories have to be derived within 
the CTP formalism and the presentation of the dynamics of the electric current. 

In Section II, we point out that inertial forces show a surprising similarity with genuine interac¬ 
tions. The equations of motion become highly involved and lead to non-trivial dispersion relations 
for the collective excitations. For us, the non-linear transformation of interest in the electron 
gas is the field transformation from that of the electron, -(/^(x) and V>(x), to the Noether current, 
j^{x) = 

In Section III, we proceed with an introduction of the CTP formalism. Special emphasis is 
placed on the structure of the CTP propagators. Physical phenomena that arise in the CTP 
formalism from the coupling between the two time axes, namely decoherence and irreversibility, 
are discussed in Section IV. We also discus why decoherence and irreversibility can be present in 
harmonic systems. 

The CTP scheme is then applied to the electron gas in Section V by working out the effec¬ 
tive action for the current at the quadratic level and by using the one-loop approximation. The 
corresponding equations of motion are presented for longitudinal and transverse components of 
the current. Section VI is devoted to the discussion of the infrared, hydrodynamical limit of the 
solutions. We first consider the linearised equations of motion and decoherence in a dense ideal 
gas with vanishing temperature. We then couple it to a heat bath and finally turn on the Coulomb 
interactions. 

Finally, we present our conclusions in Section VII. Three Appendices include a brief discussion 
of the linear response formalism, the structure of the free Schwinger-Keldysh propagators and the 
summary of the one-loop current-current Green’s function calculation. 


II. INERTIAL FORCES AND INTERACTIONS 

We begin this paper by analysing the dynamics of physical systems in which the quantities of 
interest are non-linear combinations of the fundamental degrees of freedom. 


A. Particle dynamics 


To examine non-linear field transformations, consider first a simple harmonic oscillator, 


T 1.2 1 22 

L = -mx - moj x . 

2 2 


( 3 ) 


We can introduce a non-linear coordinate y = x^/2, under which the Lagrangian becomes 


1 f m\ .2 1 1 2 2 


( 4 ) 
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where ojy = 2uj. The system in Eq. (4), expressed in terms of y, appears to have a position- 
dependent effective mass, m{y) = ml2y, a periodically diverging speed and oscillations, which are 
unable to pass through the point y = 0. Despite the fact that both systems describe the same 
physics, the inertial force exerted by the coordinate-dependent mass of the oscillator y{t) hides the 
simplicity of the motion expressed in terms of x{t). 

The lesson to be learned at this point is that inertial forces, arising from a non-linear change of 
the coordinates, may disguise the normal modes of a harmonic system beyond recognition. Let us 
suppose that our system, described by the coordinate x, obeys dynamics described by a quadratic 
action, S[x], and that we are interested in the effective theory of some non-linear combination of 
the original coordinate, y = y[x]. The action, ^[^[y]], is then non-linear and contains terms that 
make the physical system appear to behave as an interacting system [26]. 

B. Classical fields 

Let us increase the complexity of the model and consider a complex scalar field, 0(x), governed 
by a translationally invariant quadratic action, 

S = cl)*C^ + + + k^. (5) 

The products of the type fg denote space-time integrations, 

/y = ^ y d^xf{x)g{x) = c j ■^^f{-q)g{q), (6) 

with the Fourier transform defined as^ 

f{q) =^J (7) 

The source k{x) is coupled to the composite field <l>(x) = 4>*{x)4>{x) and the kernel C{idy) can 
be a polynomial in space-time derivatives. The normal modes of the model (5) with j = k = 0 
are plane waves with a wave-vector = {uJo{q),q), where C{qy) = 0. If we turn on a source 
j with four-momentum then all of the energy-momentum gets absorbed exclusively by the 

normal mode q^. The system’s response to an external source j thus reflects the simplicity of a 
harmonic system in which each normal mode remains independently excited. On the other hand, 
the energy-momentum injected through the source k, which is coupled to the composite field <I>, 
gets transmitted through the system in a completely different way. It spreads over all normal modes 
with the energy-momentum conservation being the only restriction. The dispersion relation of the 
response is modified in a highly non-trivial manner compared to the previously excited normal 
modes with coo{q). 

The above statements can also be understood in the following way: after turning on the source 
k, we not only lose the original normal modes but the response of the system begins displaying 

Throughout this work, we will be using the rjin, = diag{-|-l, —1, —1, —1} sign convention for the Minkowski metric 
tensor. 
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correlations among those modes. It is important to stress that correlations are usually thought of 
as the hallmark of genuine interactions. In this language, we can understand the non-interacting 
nature of the model diagnosed by j as the absence of correlations among the normal modes. In 
fact, the external source j{q) ~ 6{q — p), which is local in Fourier space, induces a local response, 
6{q — p). This is no longer the case when the source k is used to diagnose the dynamics 
since k{q) ~ 5{q — p) induces a response in the whole Fourier space because the momentum p^, 
injected by the source k, is spread over inhnitely many plane wave normal modes of the k ^ 0 
model. 

It is also instructive to compare the time evolution of <I>(x) induced by a wave packet in k(x) 
with the spread of the wave packet in quantum mechanics. In standard quantum mechanics, the 
spatial Fourier transform of the probability density can be written as 

p{t,P) = J + ( 8 ) 


where V'(q) denotes the Fourier transform of the wave function at t = 0, assumed to be regular and 
localised in Fourier space. The probability density is also localised at t = 0 and its p-dependence 
displays a convoluted peak. As time passes, the oscillating phase factor makes this peak more 
and more localised in Fourier space, leading to the spread of the wave packet in position space. 
The time reversal, being an anti-unitary transformation, leaves the Schrodinger equation invariant 
and the spread of the wave-packet results from the choice of the initial condition rather than the 
breakdown of time reversal invariance. It is possible to find initial conditions for the wave function 
with a singular phase in Fourier space so that the wave packet becomes narrower in time. 

The field <h(x), induced in the field theory (5) with j ^ 0 and k = 0 is given by 


f ^ {p+q)t 


ah 


(27r) 


Res C ( (q) , q ) Res C ( (p + q), p + q 




X J (q), q) j (p + q),R + q), 


(9) 


where Res (q), q) is the residue of the retarded Green’s function and the summation extends 

over the poles. The comparison with Eq. (8) indicates that if the source j has a regular and localised 
Fourier transform, then the peak in <h(t,£c) spreads as the time passes. The interference among 
the independent normal modes leads to diffusion-like processes as in the case of Landau damping 
[27]. As a result, the mixing of infinitely many normal modes may make the effective dynamics of 
the composite field, <h(x), appear to be diffusive. 


C. Quantum fields 


To continue with our chain of increasingly complex models, consider a charged scalar quantum 
field. 


d^q 


■''(*) = E y ^ [®(“o) +e(-“o) '’^(-9) 




( 10 ) 
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where the operators ac{q) and bc{q) annihilate a particle and an anti-particle, respectively, and the 
sum is over the different particle modes of the Lagrangian (5). The non-interacting nature of an 
ideal gas is clearly apparent on the level of the one-particle Green’s functions. The factorisation of 
the higher-order Green’s functions of the elementary field. 


G2n — ( 0 


r 


|(/)(xi) • • • (f){Xn)4)Hyi) ■ ■ ■ </>^(yn)} 


( 11 ) 


according to Wick’s theorem, implies that all connected Green’s functions with n > 1 vanish. In 
fact, there are n elementary excitations contributing to G 2 n, each of them controlled by a pair of 
particle or anti-particle operators, ac(q)aj(q') or bc{q)bl{q'). They can be represented by propa¬ 
gator lines in the Feynman diagrams and G 2 n is the sum over all possible different combinations 
of the pairings. Since the elementary excitations are non-interacting, their amplitudes factorise for 
all the combinations. 

The Green’s functions of the composite operator <h(x), 


^2n = (0|r{$(xi)---$(x„)}|0), 


( 12 ) 


have a much more involved structure. In fact, each composite operator is made of two elementary 
excitations and therefore <h(x) can be paired with two elementary field lines in a Feynman diagram, 
as for example in Figures 1 and 2. In other words, the measurement of ‘h(xi) generates two 
elementary excitations, which have to be removed by other operators appearing in Fl 2 n- The lines 
drawn between the composite operator insertions can therefore give rise to connected diagrams at 
arbitrary orders. As a result, the simplicity of the non-interacting, ideal gas becomes completely 
disguised by the composite operator analysis. 



FIG. 1. 


An example of a scalar two-point Green’s function ring diagram for a composite operator. 



FIG. 2. An example of a scalar n-point Green’s function ring diagram, which contains n current insertions 
of ® along the ring with propagators connecting the neighbouring pairs of the operator insertions. 

Even though both the non-interacting and the interacting system share the property that they 
support non-trivial connected composite operator Green’s functions of arbitrarily high orders, there 
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is a marked difference between these Green’s functions. For an ideal gas, they are given by a single 
Feynman loop-integral. In the presence of interactions, on the other hand, there exists an infinite 
series of distinct Feynman diagrams. For an example of a diagram with added interaction vertices, 
see Figure 3. Let us consider a composite operator that is an A^-th order homogeneous multinomial 
of an elementary free field. It is easy to see that in the non-interacting case its Green’s function 
with n legs, corresponding to an n-point function, is given by a single Feynman diagram with 
n{N — 2)/2 -|- 1 loops. For example, all n-point functions for a bilinear operator with N = 2, as 
depicted in Figure 2, are given by one-loop diagrams. The numerical values of the many-body 
correlations can therefore be as involved in a non-interacting gas as in an interacting system. 
However, the elementary processes that build up these correlations are far more complicated in the 
interacting than in the free case. This is consistent with the findings of Section VID, where we 
show that the Coulomb interaction re-summed propagator only modifies the numerical values of 
the coefficients in the linearised equation of motion for the current, compared to those derived for 
an ideal gas. We note that our analysis will neglect the vertex corrections as well as the electron 
self-energy, which would qualitatively change the behaviour of the gas. We will discuss these issues 
in detail below. 



FIG. 3. A simple two-point current-current Green’s function diagram in the presence of interactions. 


To conclude this section, let us consider a particle described by a complex field which is 

invariant under the phase transformations (/>(x) —)■ e''^cj){x) and cj)\x) —?• e~''^(j)\x). Its interactions 
can be described by vertices constructed out of the bi-local field <I>(x, y) = (l){x)cj)'^ {y). Furthermore, 
^{x,y) also generates the field’s excitations. In a sense, this is a generalisation of the ideas on 
bosonisation because the excitations are controlled by a bosonic bi-local field, both for bosons and 
fermions. The dynamics of such excitations is non-local. However, by adopting the idea of the 
operator product expansion, it can be characterised by a set of infinitely many local fields, 

— (j){x + z)(j)\x - z)|^=o- (13) 

Ozl^l * * * Ozl^n 

One expects that the long wavelength excitations of the system are generated by the fields with 
small n. It is therefore reasonable to look for a simple effective infrared theory in terms of <ho- 
Such effective dynamics may contain dissipative forces because the excitations that are controlled 
by the fields <hn, with higher re, represent an environment for the long-distance modes. Note that 
this environment is well defined even though it has no corresponding direct product structure in 
the Fock space. 
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III. CTP FORMALISM 

The standard single-time axis formalism of quantum field theory, developed for the calculation 
of transition amplitudes between pure states, is not sufficient to describe the time evolution of 
expectation values. To overcome this problem, we will use the Schwinger-Keldysh, Closed-Time- 
Path (CTP) formalism [13] in order to derive the equation of motion for the expectation value of 
the Noether current j^. This powerful formalism has been applied to numerous calculations in 
condensed matter physics as well as in high energy physics, see e.g. [28, 29]. To make the paper 
self-contained, we will devote this section and Section IV to the introduction and summary of the 
relevant CTP techniques employed in this work. 


A. Generating functional 


To make our presentation of the CTP techniques as simple as possible, we will work with a 
specihc example. We will consider the dynamics of a neutral scalar field 4>{x) during the time 
interval ti < t < tf, defined by the Hermitian Hamiltonian density H(x). In standard QFT, 
transition amplitudes between pure initial and hnal states are calculated by using the generating 
functional 


en 


w\j] ^ 




-r f d‘^x[H{x)-j{x)(j){x)] 




(14) 


The source j{x) is introduced so that we can generate n-point Green’s functions and facilitate 
the perturbation expansion. In order to compute expectation values and allow for an evolution to 
(and from) mixed states, one generalises the expression (14) to the trace of the final density matrix 
pj = p{tf), written in the Heisenberg representation as 


giW[j+,j ] ^ 


r 






^ f d'^x[H{x)+j {x)<j>{x)] 


(15) 


where pi = p{ti) stands for the density matrix of the initial state and T* denotes the anti-time 
ordering. The sources j^{x) generate observables through functional differentiation, 6/6j^{x), 
and are set to the physical value, j~^{x) = —j~{x) = j(x), at the end of the calculation. One 
can distinguish between the two time axes by introducing two time variables C*“ and t~ for the 
time-ordered and the anti-time-ordered products, respectively. Furthermore, it is convenient to 
dehne the extended time ordering, T, in such a manner that it acts as T on C*“ and T* on t ~, while 
placing each t~ after t^. The result is the condensed expression. 


^TiWlj+J ] ^ rjY 




^fd'^xlHlx )-H{x+)+j+{x+)<l){x+)+j {x )(j){x )] 




Pi 


(16) 


where {t, x) goes from (C*“, a;) to {t~ ,x). Moreover, Wick’s theorem also becomes explicitly available 
and can be used in calculations. 

The path integral representation of (15) is 


1 = J D[(f)'^]D[4> ]pi[(j)'^{ti,x),4i (tj, \ (17) 
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where S'[(^] is the action, which includes Feynman’s ie prescription. In this paper, we will make 
use of the notation 

= J (18) 

where the CTP doublets are (p = ( 0 +, 4)~) and j = and the action is 

= (19) 

S* here denotes the complex conjugate of S. The CTP symmetry, 

s[rA-^] = -s*[cp+,ri ( 20 ) 

plays an important role in restricting the structure of the Green’s functions and effective actions. 

The unitarity of the time evolution, on the level of the whole system, is expressed by the 
preservation of the total probability. This is reflected in the fact that the trace of the density 
matrix, (15), calculated for a physical source, = —j~ = j, equals to 1, or equivalently, 

W[j,-j] = 0. (21) 


Is the CTP scheme simply a trivial generalisation of the usual, single time axis formalism? 
The generating functional (15) factorises into the product of the conventional transition amplitude 
(14) and its complex conjugate, with W[j,—j] = W[j] — if the initial state, iTj), evolves 

under the influence of the given external source into f) at the final time. Such a trivial relation 
between the single and the double time axes cases is no longer valid when either the initial state 
is mixed or we use “non-physical” sources, j~^ ^ —j~. In such cases, the CTP scheme covers new, 
non-trivial physical phenomena, which could not be described by (14). The essential point is that 
in calculating the expectation values of operators, we always encounter non-physical sources, such 
as 


Tv[(j){x)pi] 


SW[j] ^ 6W[j] 


( 22 ) 


which explains the need for the CTP formalism, as for instance in the theory of linear response. 
We will make use of this feature in the construction of the effective dynamics for the expectation 
value of the electric current, described in Section V. As mentioned above, another circnmstance 
that makes the CTP formalism a necessity is the presence of a mixed initial state, which happens 
when the system under consideration is open due to its coupling to an environment. In Section 
IV, we will explain that this extension is in fact essential for effective theories. 


B. Propagator 


The CTP formalism with its doubling of the degrees of freedom allows us to set up perturbation 
expansion for retarded Green’s functions. They are completely encoded by the CTP propagator. 


iD^^'{x,y) = Tr r |(()'"(x)(y) | p* 


(23) 






13 


defined by the generating functional (18). We will use a to denote the CTP indices + and —. The 
generalised time ordering, T, agrees with the usual one, T, on the positive time axis. Therefore, 
D++ is the Feynman propagator, 


Tt [r {0+(x)</.+ (y)} |0)(0|] = (0|r{</.(x)</.(y)}|0). (24) 

The action of T is trivial if the two operators belong to different time axes, giving us 

Tr[r{r(x)<^+(2/)}|0)(0|] =Tr[0(x)</.(y)|O)(O|] = {0\cPix)cPiy)\0). (25) 

Hence, the off-diagonal components of the propagator are the Wightman functions without time 
ordering. The remaining components of D°'^'{x,y) can be found by complex conjugation, leading 
to the block matrix form 

), ( 26 ) 

Wxmy)) {Timmir) 

The CTP propagators for free bosons and fermions are given in Appendix B. The CTP identity. 


T{A{tA)B{tB)} + T* {A{tA)B{tB)} = A{tA)B{tB) + B{tB)A{tA), (27) 


valid for bosonic operators, restricts the propagator to the standard CTP form, 

(d^ + W^ -Df+ iD^\ , , 

D={ . . , 28 

\D^ -D^ + iD^J 

where the functions D'^, and appearing in the matrix elements are real. In the bosonic case, 
the exchange symmetry (u, x) -f-)- imposes D'^{x,y) = D^{y,x), D^{x,y) = —Df{y,x), and 

D\x,y) = D^{y,x). 

The Fourier transform of the Wightman function, 

zD-+{p) = eip^)S{p), (29) 


is the spectral function of excitations, which are generated by 4>{p) in a system with translational 
invariance and S{p) > 0. The relation (29) allows us to express both and in terms of the 
spectral function, which leads to the spectral condition 


D^p) = sgn{p^)iD\p), 

and the CTP propagator can thus be specified by only two real functions, 

+ ®Sn(p°)T)7(p) -20 (-p°)T>7(p) \ 

^ y 2Q{p^)D^{p) —D^{p) + sgTi{p^)D^{p)j^ 

where the positive definiteness of the norm imposes the bound 


(30) 


(31) 


iQ{p^)D^{p) > 0 . 


( 32 ) 
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The inverse of the propagator (28) is given by 


b-^ = a 


+ -A^ + iA* 


cr, 


^A-^ + iA* -A” + iA7 ’ 

where ct is a diagonal “metric tensor” of the form a = diag(l, —1). Furthermore, 

A’'’“ = 

A* = -A’'L»*A“, 


(33) 


(34) 

(35) 


where A"'(x,y) = A"’(y,x), A'^(x,y) = —A'^(y,x), A*(x,y) = A*(y,x), A'" = A"’ + A-^ and 
A“ = A” — A-^. We also note that the spectral condition, i.e. Eq. (30), yields 


A'^(p) = sgn(p°)iA*(p). 


(36) 


Even though the preceding remarks apply to interacting fields, it is instructive to consider free 
fields in a harmonic model given by the action 


S,ff[^] = - [ d^x (</.+ , 0 -) ( ^ . 

2j M-A/-fA* -A’^ + zA* 


(37) 


The external source, introduced in the generating functional (15) generates a non-trivial expecta¬ 
tion value, which can be obtained from either one of the two time axes. 


(^w) (;)=-/Ao(-.i/)T“)). 


(38) 


showing that D” = and can be identified as the retarded and the 

advanced Green’s function, respectively. Since these Green’s functions are real in position space, 
D^{p) = ^D^{p) and D^p) = f9T»^(p). 

IV. COUPLING TO THE ENVIRONMENT 


In the above section, we mentioned that the GTP scheme is important for the analysis of open 
systems. In this section, we will address the question of how precisely this formalism is able to 
describe the coupling between a system and its environment. 


A. Effective theories 

Let us suppose that the states of a closed dynamical system correspond to the linear space 
V. = Tis 0 Tde, and that we are interested in the expectation values of observables acting only on 
the system factor space, Us, assuming that the initial state is pure, iTj) G T-L. It is well known that 
the environment, represented by the factor space Tie, decouples and that the expectation values 
can be obtained within T-Lg if the state is factorisable, |4'(t)) = |</>(t)) ® IV’(t))) with \4>{t)) G T-ls, 
and IV’(t)) £ T~(-e- This scenario is certainly applicable when the initial state is factorisable and 
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there is no interaction between the system and its environment. But as soon as the system and 
the environment begin to interact in some manner, the factorisability of the state vector is lost. 
The system and the environment must be in a pure, entangled state; hence the expectation values 
within the system can no longer be reproduced by pure states in %s- In this case, the system state 
becomes mixed and must be represented by the reduced density matrix. 

The traditional framework to discuss effective dynamics is based on the single time axis for¬ 
malism and the generating functional (14), which is designed to describe the transition amplitudes 
between pure states. As a result, it is very difficult to use such a formalism to describe mixed 
states and open systems. An example of such difficulties is the necessity to perform an involved, 
partial re-summation of the perturbation series to cancel the collinear divergences in the effective 
theory of point charges [30-33]. This example will be discussed again below, in Section IV D. 

Another property of effective theory, which is not covered by the single time axis formalism is 
irreversibility. This is because the variational equations are conservative. In the CTP formalism, 
however, the doubling of the degrees of freedom allows us the realise irreversibility within variational 
equations by coupling the two time axes. We will further discuss irreversibility in Section IV E. 


B. Coupling of the time 2Lxes in effective theories 

The distinguishing feature of the CTP formalism is the coupling of the two time axes, which 
was implemented on the level of a phenomenological effective action for hydrodynamics in [9] . Such 
couplings are always present and the simplicity of the action (19) is misleading as it does not reflect 
the fact that the dynamical variables on the two time axes are correlated at t = tf. In fact, the hnal 
condition, 4>'^(tf,x) = does not fix the dynamical variables as in the single time axes 

path integral scheme. Instead, it establishes a correlation among (j)'^{tf,x) and (p~{tf,x), leaving 
the actual value undetermined. Furthermore, if (j)~^ and 4>~ are introduced as two independent 
variables, then the constraint representing the final-time condition implies a non-trivial coupling 
between the time axes at t = tf. 

The explicit appearance of the condition x) = x) in the path integral (17) violates 

translational invariance in time. To recover this symmetry, together with the diagonal structure of 
the Green’s functions in frequency space, one can take the limits of ti —>■ —oo and t/ —)• oo. One may 
expect that the coupling of the time axes, generated at t = tf, disappears in this limit. However, 
a more careful calculation shows that the time axes remain coupled by a time-translationally 
invariant, infinitesimally strong coupling [34]. Such a transmutation of the boundary conditions is 
actually the origin of the 0{e) terms in the quadratic action of a free particle, given by Eqs. (B3). 

In effective theories, the time axes are correlated by finite couplings. As an example, let us 
consider the system of two interacting fields, 4){x) and ipix), governed by the action 

S[(l),1p] = Ss[4>]+Sse[4>,'lp], (39) 


where Sge describes the interaction between the system and its environment. The CTP generating 
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functional for the effective theory of <^(x) is 

= J D[(f>]D[il;]expi^^S[(j)'^,'il;^] - ^S*[4i~ , (40) 

where 4>'^{tf,x) = 4>~{tf,x) and ijj~^{tf,x) = 'ijj~{tf,x). After integrating out the fields the 
bare effective action can be written in the form of Eq. (18) as 

= Ss[cP+] - 5:[<()-] + SinM, (41) 

where the last term on the right-hand side, the influence functional [35], contains the effective 
interactions generated by the environment. It is dehned by the expression 

^ J . (42) 

The form (41) of the action does not separate terms that play different roles in the effective 
dynamics. To distinguish different contributions to the physical properties of the system, it is more 
illuminating to write 

Sinfi m = Si [0+] - SI [r] + S2 m , (43) 

where ^ 2 [</>“'“, 0] = 52[0, 4>~] = 0. Si is then the self interaction and S 2 the entanglement functional. 
They contain the contributions from the single and the double time axes, respectively. Note that 
the symmetries that act on both components of the CTP doublet identically and simultaneously 
are preserved by Seff, Sinfi and 52. 


C. Perturbative self-energy and the entanglement functional 


The physical origin of the separation of the effective vertices of the influence functional into 5i 
and S 2 is easiest to see in the perturbative construction of the effective dynamics. Let us write the 
action for i/j as 


5se [</>,'!/>] = V + 5se,mt[0, V']- 

The perturbative expansion of the influence functional is then defined by 


(44) 




1 5 


’ i 6J+ 


-s: 


1 5 


i 6J 


exp<j -Wo[J] 


(45) 


where 


eiWo[J] ^ ^ 


(46) 


is the free generating functional for the held if. 

The couplings between the time axes arise from the free CTP propagator. Do, given by Eq. 
(26), which has non-vanishing off-diagonal blocks. The propagator lines in the Feynman diagrams, 
representing the perturbation series (45), that connect elementary vertices on different time axes 
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are the ones encoded in these ofF-diagonal blocks. To find the physical origin of the off-diagonal 
propagator components, write Eq. (25) as 

T)"+(x, y) = ^ (n|(/>(x)(?i(y)|0) (0|n) = ^ (0|(/>(x)|n) {n\(j){y)\0), (47) 

n n 

where {|n)} is a complete set of basis vectors. D ^ thus receives contributions from one-particle 
states, which together contribute to the trace in the generating functional (14) and the propagator 
(23) at t = tf. The internal lines of the Feynman diagrams that connect the time axes therefore 
correspond to the excitations of the field at the final time. 

The states that contribute to the trace and realise the couplings between the time axes are on the 
mass-shell because the Wightman function in Eq. (25) contains no time ordering. S 2 is therefore 
suppressed in the limit of —)• 00 if the excitations of the initial state, involving '(/^-particles, have 
a gap. 

The single time axis contributions to the effective action, ^i, describe the type of system- 
environment interactions that can be successfully handled by the conventional, single time axis 
formalism. Hence, the separation (43) splits the system-environment interactions into two parts, 
one that keeps a pure state pure and another that generates entanglement. 

D. Conservation laws 

The separation (43) of the effective interactions is particularly clear in classical systems where 
the CTP action has an infinitesimal imaginary part. The full theory usually possesses some contin¬ 
uous symmetries, which are realised in the same manner on each of the two time axes, for instance 
external space-time symmetries or internal global symmetries. Effective vertices in Si are called 
conservative because they correspond to a traditional, single time axis action and thereby preserve 
the conservation laws for all conserved currents of the full theory, derived from the Noether the¬ 
orem. The functional S 2 , which mixes the two time axes, contains the remaining, non-conserving 
and dissipative system-environment interactions, which make the system open [36]. 

Such a formal separation of the CTP effective vertices, based on conservation laws, is less clear 
in the quantum case where the effective action has an imaginary part and '^Si and 352 also 
contribute to dissipative phenomena. The system-environment interactions generate two kinds 
of dissipative phenomena. Firstly, the finite life-time of quasi-particles due to the leakage of the 
pure quasi-particle states into the environment, called damping, is expressed by Q'Si. Secondly, the 
system-environment entanglement makes the quasi-particle states in the effective theory mixed and 
generates decoherence, which will be discussed in Section IV F. The latter phenomena are encoded 
by 3'5'2. 

In the effective theories of charges in QED, the typical example of the effects described by Si 
is vacuum polarisation, which is produced by the virtual soft photon content of charged particles. 
Virtual intermediate excitations cannot violate conservation laws in the asymptotic, on-shell states. 
The finite energy resolution of observations, performed in a finite amount of time always leaves the 
on-shell, soft photon content of charged particle states unresolved and the realistic charged states, 
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which avoid the collinear divergences of non-degenerate perturbation expansion, are defined with 
a small but finite energy spread. A typical effective interaction described by S 2 is therefore due to 
real soft photon contributions to charged particle states. 

The cancellation of the collinear divergences of the non-degenerate perturbation expansion in 
the effective theory of point charges, mentioned above, requires us to sum up the real and the 
virtual soft photon contributions. This is a very natural procedure within the CTP formalism, 
as the former and the latter are automatically encoded by Si and S 2 , respectively. The same 
construction in terms of the single time axis formalism is rather involved but is required to compute 
finite transition amplitudes in the effective theory. 

It is essential to note that the mixing of fields in S 2 breaks space-time translational symme¬ 
tries, thereby causing the physical system on each of the time axes to break energy-momentum 
conservation. This is the way an effective action that encodes the dynamics of only a subset of all 
degrees of freedom can describe dissipation, i.e. loss of energy from the relevant effective degrees 
of freedom to those that were integrated out. This fact was used for example in [9] to write down 
a classical dissipative action of a fluid with bulk viscosity. 


E. Irreversibility 

It is important to distinguish between the loss of time reversal invariance and the irreversibility 
of an effective theory [37]. By assuming that the microscopic dynamics is time reversal invariant, 
the only source of time reversal non-invariant effective interaction is the initial condition of the 
environment. The time inversion odd piece of the action (37) is A-^, which completely determines 
S 2 according to the spectral condition (36). Therefore the harmonic theory (37) breaks time 
reversal invariance if ^ 0. The form (B3) of the inverse free scalar propagator shows that 
an infinitesimally strong violation of time reversal invariance is sufficient to produce finite time 
reversal odd expectation values, e.g. / 0 in Eq. (38). 

The amplification of infinitesimally small symmetry breaking, on the level of the action, to a 
finite symmetry breaking by expectation values is the hallmark of spontaneous symmetry breaking. 
Let us again consider a free field, i.e. in the model (37) with infinitesimal A-f, given by (B3). Such 
an amplification takes place for a free field because of the inversion D = A“^, with its sole role in 
the dynamics being the representation of the initial conditions. This is no longer the case if A^ 
assumes a finite value. S 2 is then finite and the dynamics, governed by the equation of motion 
A(f> = —j, is open and non-conservative. An environment with a sufficiently large capacity for 
absorbing energy and a non-trivial spectral weight at vanishing energy then makes the effective 
dynamics irreversible [38]. 

The microscopic realisation of irreversibility is dissipation, which is the leakage of the system 
into the environment. The resulting effect, i.e. damping, is described by QSi and 5R52, with the 
former being responsible for the finite life-time of the excitations and the latter contributing to 
diffusive forces in the equation of motion derived from the effective theory. In a classical effective 
theory, such as the one considered in [9] , the only way to introduce dissipation is through S 2 ■ This 
is because the imaginary part of the Lagrangian is infinitesimally small. 
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The emergence of a complex action in a conventional, single time axis path integral signals 
non-unitary time evolution. However, the unitarity condition, expressed by Eq. (21), is preserved 
in the effective theories, which therefore always preserve the unitarity of the microscopic theory. 
Such a global unitarity is maintained by compensating for the leakage of the system state into the 
environment, described by with the system-environment entanglement expressed by 82 - 


F. Decoherence 


Decoherence [39, 40] is the suppression of the off-diagonal density matrix elements and is a 
basis-dependent phenomenon. In fact, the density matrix, being Hermitian, is always diagonal in 
a suitable basis. Here, we wish to consider the decoherence of a local operator 4>(a;) = F[(f)(x)]. 
More precisely, we are interested in the suppression of the off-diagonal elements of the reduced 
density matrix for <I>(a;) in the field-diagonal basis where the basis vectors are the eigenvectors of 
$(a;). It is convenient to use the Fourier transform, ‘h(fc), rather than the coordinate-dependent 
field because the eigenstates of this operator contribute independently to the suppression in our 
truncated effective theory. 

For the purposes of analysing decoherence, it is useful to generalise the CTP formalism, which is 
aimed at computing the trace of the density matrix, to the Open-Time-Path (OTP) scheme, which 
is able to reproduce the full density matrix. Instead of Eq. (15), one defines the OTP generating 
functional as a matrix element of the density matrix, identified by the spatial configurations <p^{x) 
in the presence of external sources. 


dW[jFr;>f>r^<f>7] = o\ 6 +. 6 - 


)4’ j'tf] 




-j:fd-^x[H{x)-j+ix)(j){x)] 


}piu)r*{ 


jrjFx[H{x)+j {x)(l>(x)] 


e'i 


'/’/■)) (48) 


with finite tf. The path integral expression for this functional is 






^f{x)=4>(tf,x) 


where the action satisfies the OTP symmetry, 


(49) 


= -S*[4> ,0^,</>/,</>/]■ 


(50) 


The generating functional for the reduced density matrix of a composite operator F[(p(x)], the 
OTP analog of the effective theory, is then 

(51) 

where the state |4>(a;)) is an eigenstate of F[4>(x)] with an eigenvalue <h(a;). In the path integral 
language, the expression (51) can be written as 


efi 


kmjAf] = 


j T»[(^]efi 




ll6Uf{x)-F[4>^{tf,x)]), 


(52) 
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or in terms of an effective action, 

I D[^e^nS.ffm+m4'], (53) 

f ,x)=^ f {x) 

The effective action in Eq. (53) is given by 




—ihlog 


X 


(54) 


which has the structure of Eq. (43). If F[(j){x)] does not determine (j){x) uniquely or if the initial 
state, described by p{ti), is mixed, then the path integral on the right-hand side of Eq. (54) 
becomes non-trivial and Sinfi ^ 0. 

Decoherence of F\(j){x)\^ i.e. the suppression of the off-diagonal reduced density matrix elements 
for the composite operator, is governed in the OTP formalism by 4*“]. It is a property 

of the state at each given instant, which gradually builds up through temporal evolution. Its 
dynamical origin is consistency [41]. This is the suppression of the contributions to the density 
matrix that correspond to well-separated OTP doublet trajectories, where <I)'''(x) — ‘h“(x) assumes 
significant values in large regions of space-time. Such a build-up of the suppression can be detected 
without having access to the full density matrix. In fact, decoherence must also appear as the 
suppression of the CTP doublet trajectories that are well separated for a significant amount of 
time. In our calculation, we will give up the details provided by the access to the off-diagonal 
density matrix elements of the OTP formalism. Instead, we will use the simpler CTP effective 
theory and only compute the trace of the density matrix to analyse decoherence in the electron 
gas. 


G. Generalised fluctuation-dissipation theorem 

Quasi-particles are dressed by the conservative interactions and have a finite life-time if 
makes the effective theory irreversible, as explained in Section IV E. Non-conservative interactions 
make the final state mixed even if the initial state is pure and 3'S'2 generates decoherence, cf. 
Section IV E. A common element of these phenomena is that both irreversibility and decoherence 
result from the quasi-degeneracy of a continuous spectrum. Such a common origin is actually 
rooted in a more fundamental identity of the CTP two-point functions, which can be considered 
as a generalisation of the fluctuation-dissipation theorem. 

The fluctuation-dissipation relation is usually discussed in the context of a linear equation of 
motion. Hence, we consider a translationally invariant, harmonic effective action given by Eq. (37) 
where the quasi-particles are defined by the single time axis action 

‘S'oi'/’j = ^ j (A" iA*) 4). (55) 

The entanglement with the environment can described by the correlation functional 

Scorr[^] = f d'^Xcjfl' ^A'^ — zA*^ (/)“. 


(56) 




21 


We will use the Keldysh basis [17], = (/> zb in which 


Seff[^] 



d}x 




(57) 


We also introduce the parametrisation = j/2 ± of the external sources, where j^{x) is the 
physical source. This represents the experimental setup in which the system is driven adiabatically 
from the perturbative vacuum at ti = —oo to a desired initial state at t = 0. The source term is 
thus 


j^ = j(l) + jP(j)'^, 


(58) 


indicating that the expectation value {4>{x)) is generated by the book-keeping source j{x), which 
should be set to zero after the functional derivatives had been taken. 

The translationally invariant operators commute and each CTP block is diagonal in Fourier 
space, a property which simplifies the discussion. The common origin of irreversibility and deco¬ 
herence in this simple model is the fact that the eigenvalue A*(g) plays two different roles in the 
effective dynamics. On the one hand, it is the imaginary part of the inverse Feynman propagator 
of the action (55) and as such it describes the inverse life-time of the plane wave q, 

A\q) = l/T{q). (59) 


It serves as a measure of strength of the breakdown of time reversal invariance. Note that A^{q) > 0 
as a result of the positive definiteness of the norm in the effective theory, cf. Eq. (32). On the 
other hand, A®(g) enters into the imaginary part of the double time axis action and 


eS'S'aM 


exp 


1 f dS 

2J (27r)4 


A^iq)\AQ)f], 


(60) 


indicating that the off-diagonal matrix elements of the density matrix are exponentially suppressed 
in (j)'^ = (j)'^ — (f)~^ and that the width of the corresponding Gaussian is 1/A*(g'). In other words, 
T{q) is a measure of the contribution of the plane waves to decoherence. It is worthwhile noting 
that a lesson of Eq. (36) is that both decoherence and time reversal invariance arise when the 
time-translationally invariant couplings of the two time axes have a finite strength. 

The coupling between the time axes, A-^, may be local in time. However, the decoherence of 
the state builds up in time in a non-local manner. In fact, sgn(a;), which appears in Eq. (36), 
contains slow time dependence. 


—e sgn(w) = -P —, 


27r 


Tit 


(61) 


if the integration is carried out over t. We use P to denote the principal value of the integral. 

What happens in an interacting system? The form (31) of the two-point function, which is 
also valid for interacting, composite operators, shows that the imaginary parts of the diagonal and 
off-diagonal blocks are given by the same real function, i.e. the far field component. The structure 
of the action (37) of our simple harmonic model therefore also remains present in the quadratic 
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part of any effective action for interacting theories. As a result, irreversibility and decoherence 
have a common origin as far as small fluctuations around a stable state are concerned. 

The double role of A*, discussed above, can be considered as an extension of the fluctuation- 
dissipation theorem. On the one hand, A*(p) is the imaginary part of the inverse Feynman propa¬ 
gator of the action (55). Furthermore, it also provides the imaginary part to the inverse retarded 
Green’s function, A-^ (p) = sgn(p^)A*(p), and thus controls the “leakage” of the quasi-particle 
states towards the environment, which is an elementary dissipative process. On the other hand, 
A““’“(p) = 20(p*^)A*(p) is the quadratic form of S 2 and therefore represents the interactions with 
the environment, which are the origin of fluctuations. Such a relation between fluctuations and dis¬ 
sipation holds even in the absence of thermodynamical reservoirs so long as the linearised equation 
of motion is reliable, i.e. the system is in a state that is resilient against microscopic fluctuations. 


V. THE ELECTRIC CURRENT IN AN ELECTRON GAS 

We are now ready to begin analysing the dynamical properties of the electric current in an 
electron gas by using the full CTP machinery introduced above. The relevant theory of the mi¬ 
croscopic physics of electrons and photons is quantum electrodynamics (QED). Our goal will be 
the derivation of the effective action for the Noether current in the quadratic and one-loop 
approximations of QED at finite density and temperature. 


A. Effective theory for the current 

We begin by introducing the microscopic theory of QED, which we will use to compute the 
Green’s functions for the current, sourced by a^, as well as the electromagnetic field, sourced by 
j^. The relevant generating functional is 




J D['ip]D['tp]D[A]exp J d'^x {P ^ A + jA |. 


The free photon propagator is 


(62) 


d^Ap) = -d{p-o)t^^^, 

with — 9^9^/n, and the free electron propagator is 


(63) 


F{p) = (^+ 


D{p; m) + 27ri(5 ( p^ — 


m?c^ 

~W~ 


n{p) 


1 - 1 ' 


(64) 


where 


D{p-,rrP) = 




-2TTi5 {f - 0(-pO 


-2T,i5 {f - 0(pO) 




(65) 
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The occupation number density n{p) is given by the expression 

Q(/) ^ 0(-p°) 

g/ 3 (ep-/x) g/ 3 (ep+Aj) _j_ I ’ 

in terms of the single particle energy, Cp = c\/m‘^c^ + and the chemical potential p. The UV 
divergences are the same as in the single time axis case, but special care is required to ensure that 
the composite operator Green’s functions are finite. The similarity between the current-current 
two-point function in the case of QED with non-interacting electrons and the one-loop photon 
polarisation tensor in the standard interacting QED implies that we require an 0(e*^) counter-term 
even for the ideal gas [42]. 

The functional VE[a,j] can be used to construct the effective action through a functional Leg¬ 
endre transform. 




W[a,j] = r[J,A]+aJ + jA, 


(67) 


where the new variables 

6 a 

are the expectation values, 


^^_ 6W[a,j] 

Sj 


( 68 ) 


J^{x) = {j^{x)) = TT[pip{x)-f>^'il;{x)], Ap{x) = Tv[pAp{x)], (69) 

for a = j = 0. 

In what is to follow, we will make use of the parametrisation 

a=^ = |±a, (70) 

for the sources. The function a represents a source that can be adjusted by a suitably chosen 
physical environment and leads to unitary time evolution, cf. Eq. (21). The other component, a, 
is only used as a formal book-keeping device. The Legendre transform of a real, convex function 
can be defined either geometrically or algebraically. We follow the latter route and use Eqs. (67) 
and (68) to define the effective action, a complex functional, in an algebraic manner. The inverse 
Legendre transform is based on the variables 

<. = -^ 1441 . pi) 

6 J 6A 

Hence, the variational equations of the effective action are satisfied by the expectation values, 
which are obtained by setting the external sources to zero. 

We are interested in the effective current dynamics and the corresponding effective action, r[J], 
which can be obtained by eliminating A from r[J, A] by making use of the second equation in (71). 
The effective action is real in the physical case, j~^ = —j~, 0 "*“ = —a~; therefore it is sufficient to 
retain only the real part, JRT, of the equation of motion. 

We write the external source coupling as 

aJ = aJ aJ^^ 


(72) 
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which yields 


SW 

6 a 


+ J ), 


da 


(73) 


showing that J is the current expectation value and = 0 in the physical case with a = 0. 

To find the linearised equation of motion, we need the quadratic parts of W and T, which will 
be calculated by expanding in h and keeping the first two orders. We start by integrating out the 
electron field, 


^kW[a3\ ^ J I)[i]exp|y d^x 
and keeping the O {h^d^) terms, giving us 




(74) 




= / D[A] exp < / d^x 


--^dGd+\kA +-^AD ^A 

In h 2th 


(75) 


The inverse of D, 

p2 

D-^ = - ^dGa, (76) 

and the source, 

k = j -\ — dGd, (77) 

c 

are given in terms of the one-loop current-current Green’s function, represented diagrammatically 
in Figure 4, 


G‘P,Zixi,X2) = -ihtr:[F^^^^Axl,X2h^,2P‘'''''Hx2,Xl)'yf,^]. (78) 

The presence of a neutralising, homogeneous background charge was assumed in Eq. (77) to avoid 
divergent tadpole contributions. The Gaussian integral (75) yields 

W[aJ] = -^dGd - ^jDoadd - ^Dq + ^DodGdDo^ j, (79) 

at the quadratic order in the sources. 

Although the orders of the expansion of VF[a,j] in h and in the number of loops correspond 
to each other in standard QFT calculations, this is no longer the case when the effective action 
is considered. The reason is that the variables of the effective action have different orders in h, 
J = O (h) and A = O {kP), as opposed to the variables of W, d = O {dP), j = O (dP)- The 
Legendre transform (67) therefore gives 

T[J,A] = ^JG-^J+^AD3 ^A-eAdJ (80) 

for the next-to-leading order, O (d) effective action. 

The Maxwell’s equation, 


A = eDQaJ, 


( 81 ) 






FIG. 4. Current-current two-point Green’s function diagram in QED. 



FIG. 5. Current-current two-point Green’s function with the Schwinger-Dyson re-summed propagator, 
can be used to eliminate the photon field so that we obtain an effective action, 

r[J] = ( 82 ) 

with the Schwinger-Dyson re-summed propagator, represented in Figure 5, giving us 

- 

C = G~^ - ^d-Doa. (83) 

The calculation of G is summarised in Appendix C. 


B. Equation of motion 


All bosonic two-point functions, including G, have the same CTP matrix structure as Eq. (28). 
This allows us to define the retarded and advanced parts of the linearised equation of motion for 
the current. 


r,a _ /^r,a\ —1 


^r,a ^ 


- — 


(84) 


where Dq’“ are given by Eqs. (C41) and (C6), respectively. The effective action, expressed in terms 
of the near and the far components of £, £"■ = (£^ -|- £“) and = (£'' — T“), is 


1 -iiC^-Cff+C- z(£--£/) 2 + £/\ /j+ 

^ J M i(£- - £/)2 - _i(£n_^/)2_^n ) I j- 


(85) 


cf. Eq. (33), which results in 

r[J, J^] = -1 ^ ^ (86) 

when expressed in terms of the Keldysh basis, i.e. Eq. (73). The form of (85) shows that ^82 
and ^82 are given by the kernels (£”■ — £7)2^ respectively. The form (86) of the effective 
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action then reveals that Q'F governs decoherence. On the other hand, the equation of motion for 
the physical expectation value J, for which a = 0, i.e. the variational equation for J'^ at = 0, 

a = J, (87) 

arises from 3f?r. We note that an effective action with the structure of this type, i.e. mixing J"*" 
and J ~, was used in [9] to describe a classical dissipative fluid. 

The Legendre transform of the effective action (82), 

W[a] = —-aCa, (88) 

is the generator functional of the connected Green’s functions for the current after the resummation 
of the ring diagrams with photon line insertions. The corresponding, improved Kubo formula for 
the induced current. 




(89) 


cf. Eq. (A14), corresponds to the inversion of the equation of motion (87). In other words, 
the equation of motion, used in this paper, could have been obtained from the traditional Kubo 
formulae for the electric current after the resummation of the ring diagrams and an inversion, 
which expresses the external source as the function of the induced response. For a review of the 
linear response formalism and its connection to the CTP formalism, see Appendix A. 

The non-relativistic limit, c —)> oo, has been well understood in QED, where the minimal 
coupling is given by the term —ej^A^/c in the Lagrangian. The source term that generates the 
electric current, is not suppressed by c. The higher-than-quadratic order terms in W, in 

the source a^, therefore represent correlations, which come from relativistic effects in QED. In 
other words, the non-relativistic mechanical equation of motion for the electric current of the non¬ 
interacting Dirac sea can contain relativistic terms in the presence of an external electromagnetic 
field. At the leading order, i.e. the quadratic level of the generating functional, studied in this 
work, there are no such unusual mixing terms. 

We can now use the non-relativistic retarded Green’s functions (G41) and the three-dimensional 
parametrisations = {4>/c,a) and = {pc,j) to find the equation of motion in Fourier space 
for the mode = {uj/c,q), 



1 

c 




+ cCt 




(90) 


where the overall factor of 1/c on the right-hand side is due to the definition (6) of the scalar 
product in Fourier space and ^ = Lo/c\q\. We are also using n = q/\q\, L = n®n and T = 1 — L, 
or in terms of the index notation, = n^nK We decomposed the above expression in terms of 
the longitudinal, 


f = 


1 





(91) 
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and the transverse part, 

= ^tT. (92) 

Finally, we are in position to derive the desired non-relativistic equations of motion, 

4> = -Lip, a = Ctj. (93) 

It is worth noting that the second equation indicates that the external source in a non-relativistic 
gas should be in the Coulomb gauge, V • a = 0. 


VI. HYDRODYNAMICAL LIMIT 

The question we wish to address in this section is whether local equilibrium can be formed in 
the presence of weak, slowly changing external sources, in space-time, leading to hydrodynamic 
behaviour of the electron gas at non-zero density, even in the absence of interactions. Somewhat 
counterintuitively, we wish to argue that this is indeed possible. While local equilibrium is assumed 
in the usual phenomenological approach to hydrodynamics, it is our goal to understand its emer¬ 
gence through a detailed derivation of an effective theory. We hope that such calculations may 
reveal new microscopic insights into hydrodynamics. To define our terminology more precisely, we 
consider the system to be in local equilibrium when the following two conditions are satisfied: 

(I) The equations of motion are local in time and contain a finitely ranged smearing in space. 
While the equations of motion of an interacting system may be local in space-time, such a 
smearing is necessary in an ideal gas due to the absence of elementary relaxation processes. 

(II) The retarded Green’s functions, describing the solution of the linearised equation of motion, 
display damped time dependence, leading to relaxation. 

In order to satisfy the condition (I), we have to restrict our attention to phenomena with 
sufficiently slow dependence on the space-time coordinates. To ensure local response, the external 
perturbations should therefore remain slow compared to the characteristic velocity. However, 
determining the precise temporal resolution for the “hydrodynamical” processes to be observable 
can be rather complicated in an ideal gas with no interactions. In the traditional phenomenological 
approach, developed for interacting systems, one assumes the analyticity of the dynamics in the 
infrared region, i.e. in a small vicinity of the point cu = 0, q = 0 in Fourier space. This renders 
the determination of the allowed space-time resolution in a hydrodynamic regime straightforward. 
Although the conservation laws tend to generate slow, long-range modes, the interactions generate a 
finite life-time, tq, for quasi-particles. This life-time acts as an infrared cutoff for the frequency and 
makes the dynamics local for w <C I/tq and |q| <C l/r^/p, where vp = hkF/m and Vmfp = vfTq 
denote the Fermi velocity and the mean free path. The absence of elementary (microscopic) 
relaxation processes in an ideal gas, i.e. as tq —)• oo, forces us to look for other sources of smearing 
mechanisms to establish local dynamics. 
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It is shown in Appendix C that the response to the external field can be conveniently expressed 
in onr case in terms of two dimensionless variables, 


X = 


muj 

h\q\kF 


and 


y = 


kp 


(94) 


The singularities of the linearised, one-loop dynamics of the ideal gas then arise at the threshold 
of the particle-hole excitations, i.e. at the border of the domain, 


T^hydr = {{y,x)\\\x\-y/2\<l] , (95) 

on the (y, x) plane. We claim that in order for the current to display hydrodynamic-like behaviour, 
uj and q must be restricted to V^ydr, where the functions (C20) and (C23) are analytic. More 
precisely, the dynamics will be most closely analogous to the usual hydrodynamics when y <C 1 
and \x\ < 1, having used Eq. (95). The resulting Laurent series for (j^) in the magnitude of 
the wave vector, |q|, starts with a non-zero, negative power of |q|. The situation here is therefore 
significantly different compared to the traditional phenomenological expansion, where {j^) is given 
by a Taylor series in a; and q. This explains the appearance of the factor l/|q| in the equations of 
motion, announced in the Introduction. 

Such a hydrodynamical regime is more restricted than that of the interacting systems, described 
by phenomenological hydrodynamics, where the analyticity of to and q allows us to consider the 
two limits w —7- 0 and q —?■ 0, independently. In the absence of interactions, the IR limit depends 
on the ratio of w and \q\, as the point to = q = 0 is approached on the (y, x) plane. To understand 
the physical meaning of the restriction imposed by V^ydr, note that x = Vph/vp, with Vph = u/\q\, 
is the ratio of the phase velocity of the external probe to the Fermi velocity. As argued above, the 
external perturbations should remain slow compared to the characteristic velocity, thus giving us a 
restriction on the size of x. As soon as interactions are turned on, the quasi-particles acquire finite 
life-time, which then acts as an IR cutoff and begins screening the l/|q| singularities. More precisely, 
it is the imaginary part of the self-energy at the Fermi surface, correcting the loop integrals (C13), 
that makes the Green’s functions analytic in q^. The traditional phenomenological approaches 
for interacting systems, which rely on the analyticity in q^, therefore correctly assume that the 
hydrodynamical regime extends to arbitrarily large values of |x|. 

We will see that the above considerations also make condition (II) satisfied for the discussed 
ranges of x and y. Details will be presented in the sections below. 


A. Spectral function at zero temperature 

The simplest insight into the response of the dynamics to an external perturbation is provided 
by the spectral weight of the excitations, given in our case by 

iG“+(g) = (0 \jy^{-q)ju{q)\ 0) . (96) 

This correlator can be recovered from the imaginary part of the retarded Green’s function in 
momentum space, i.e. the far held component. A more detailed view of the response is provided 
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by the identification of the collective modes. The dispersion relation of the collective modes, uj{q), 
of a harmonic system is easiest to dehne by the use of the residue theorem, if the inverse retarded 
Green’s function is an analytic function on the complex frequency plane. The dispersion relations 
are then defined by the zeros of the inverse retarded Green’s function, [G'’(w,q)]“^ = 0, and the 
real and the imaginary parts of u}o{q) give the frequency and the inverse life-time, respectively. A 
simplification of this prescription occurs when the frequency dependence of [G'’(a;, is linear 
and |9‘w(q)| <C |51?a;(q)|. The approximate normal mode dispersion relation is then given by 
the vanishing of q)]“^ = 0 and the inverse life-time can be approximated by evaluating 

q)]“^ on the dispersion relation. 

It is important to note that our spectral function is not analytic on the entire complex frequency 
plane at T = 0, cf. Eq. (G22), hence the residue theorem-based arguments, mentioned above, do 
not apply. One can presumably recover analyticity at non-zero temperature, however, the non- 
analytic part of the Green’s function is still present in the low-temperature expansion. The picture 
becomes clearer if we restrict our attention to a sufficiently small region near x = y = 0, where 
we do indeed recover analiticity and are able to find well-defined collective modes. This feature 
explains the universal importance of collective modes for the hydrodynamical description, as they 
can be introduced in the IR even if they become ill-defined at shorter time and length scales. 

The abundance of particle-hole states can be estimated by the spectral weight (96), given by the 
H— components of the propagators in Eq. (G3) and presented in a more detailed way in Eq. (G22). 
We depict the longitudinal and transverse spectral functions, and —Q'lS’* , respectively, in 

Eigs. 6 and 7. The numerical results are presented at metal density with kp = 10®cm“^, in units 
of Ac = h/mc = 1. The spectral weight at y ~ 0 is proportional to x and drops rapidly to zero for 
|x| = 1 — O [y). Eor y S> 1, the Eermi sphere is negligible and the holes carry negligible energy- 
momentum. Hence, the particle-hole excitations, both in the longitudinal and in the transverse 
sectors, obey approximately the same dispersion relation of a free particle and spread over the 
interval 2\x\ — 1 < y < 2\x\ + 1, for |x| ;§> 1. 

However, the spectral weights in the longitudinal and the transverse sectors are markedly dif¬ 
ferent. The majority of the longitudinal modes are present at long wavelengths and their number 
diminishes at shorter wavelengths. On the other hand, the transverse excitations are more common 
at short wavelengths. Of course, these are relative statements, as the longitudinal spectral weight 
is several orders of magnitude larger than the transverse one in the kinematical regime of interest. 


B. Excitations at zero temperature 

We begin a more detailed analysis of the current dynamics by first considering an ideal gas with 
e = 0 at zero temperature, T = 0. The retarded Green’s function and its inverse can be presented 
in the hydrodynamical limit, discussed above, by the Laurent expansion in Eourier space. The 
longitudinal retarded Green’s function, 


kpm 


( 97 ) 
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-Im T^- 


-Im S*- 



FIG. 6. Plots of the longitudinal and transverse spectral functions, shown over the (y, —x) plane, in units 
of TO = 1. (a): The longitudinal component, — {x,y), with values between 0 and 10“^. (b): The 
transverse component, —95^ {x,y), with values between 0 and 4 x 10“®. 

-Im T*~ -Im S*" 




FIG. 7. Plots of the longitudinal and transverse spectral functions, shown as functions of —oj/kp, in units 
of TO = 1. The curves with different colours represent spectral functions at different dimensionless momenta 
y = \q\/kp. In both plots the colours range from black (y = I) to progressively lighter blue colours (up to 
y = 9) in steps of 1. (a): The longitudinal component, —(w, |q|). (b): The transverse component, 
-S55+-(a;, |q|). 


is given by the dimensionless function G^, which assumes the following form in the IR limit, 

OO 

Ge= (98) 

j,k=0 


The hrst few coefficients of the series are given in Table I, with the actual order of the truncation 
of the series to be justified later, cf. the discussion after Eq. (118). The inverse Green’s function 
can be written as 




kpiTi ^ 


( 99 ) 
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y 

FIG. 8. The blue solid line shows the full, numerically computed longitudinal collective mode dispersion 
relation, X(,{q), for an ideal gas at metal density, as the function of \q\/kp. The black contour lines correspond 
to which reaches its minimum around |<7| ^ O.lkp- The plot uses the dimensionless {y,x) variables 

for the axes. 


and contains the dimensionless function 

OO 

Hi=Yl bi,j,k,o{ixyy^'^- ( 100 ) 

j,k=0 

The expansion coefficients are listed in Table II. 

As discussed in Section VI A, the real part of uj{q) at which the inverse Green’s function vanishes, 
gives the dispersion relation of the collective mode. We therefore need to look for = 0. What 
we find is a strongly damped sound wave, propagating with velocity 


= lim 
(7-s-O 


\q\ 



0.83 vp. 


( 101 ) 


Fig. 8 shows the exact longitudinal dispersion relation, Xi{q) = mL0i{q)/\q\kF, which was found 
numerically by plotting the line ^x[GW~^ = 0 for the full retarded Green’s function, given by Eqs. 
(C19), (C21) and (G41). The analytical form of Eqs. (98) reproduces the inverse of the full Green’s 
function very accurately, given by Eqs. (C19) and (C22), along with the sound wave dispersion 
relation for y <C 1. Note that the sound wave is within the hydrodynamical regime due to the fact 
that 6 £,o,o,o/^£, 2 ,o,o < 1- 

One could naively expect that the transverse sector should display no collective modes since the 
solutions of the Schrodinger’s equation only describe longitudinal waves. However, the collective 
mode spectrum, the curve 'Rx[Gl]~^ = 0, obtained numerically from Eqs. (G19), (C21) and 
(C41) and plotted in Fig. 9, suggests the presence of quasi-particles with an approximatively 
free dispersion relation, in agreement with the well known transverse waves in dissipative fluids 
[1]. This phenomenon may take place in an ideal gas by combining particle and hole plane wave 
modes, which are one-by-one longitudinal, but have a transverse component with respect to their 
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y 

FIG. 9. The plot depicts the same information as Fig. 8, except for the transverse collective modes. The 
imaginary part of the equation of motion, approaches 0 from below as |< 7 | increases. 


total momentum. For the transverse Green’s function, we find 


G 


r 

t 


Tr'^h? 


G 


t: 


( 102 ) 


with 


OO 

Gt= at,j,k,Q{ixy, (103) 

j,k=0 

in the infrared regime, cf. Table I. The absence of the term O (x*^y^) signals that there are 
transverse collective modes in the infrared regime, but that their equation of motion is non-local. 
The dispersion relations in the longitudinal and the transverse sectors are similar for y ^ 1, which 
is in agreement with the remarks made in Section VIA. 

To conclude this sub-section, we note that the plan to derive hydrodynamical equations has 
surprisingly already been partially realised in the case of a zero-temperature electron gas at finite 
density. Indeed, the longitudinal sector possesses a hydrodynamical limit, it satishes a local equa¬ 
tion of motion in the IR and displays a collective mode for |q| —?■ 0 with a linear dispersion relation, 
a; ~ \q\. The collective excitation is therefore a sound wave. However, because of the fact that it is 
composed out of non-interacting constituents, the mode is neither zero nor hrst sound. Instead, its 
existence stems from the composite nature of the operator generating it, i.e. the Noether current. 
This collective mode can thus be thought of as a composite sound mode consisting of infinitely 
many simple plane waves. To establish its existence, no additional thermodynamical considera¬ 
tions were needed as the assumption of local equilibrium is replaced by the infrared conditions 
needed to derive the local equations of motion. The transverse sector of the current has collective 
modes that can be approximated by the free dispersion relation, but their equation of motion is 
completely non-local. 
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C. Excitations at small, non-zero temperature 


In order to recover local equation of motion for the transverse modes we bring the ideal gas into 
contact with a heat bath. Because the loop-integrals in Eqs. (C13) cannot be computed in closed 
form for arbitrarily large temperatures, we will restrict our attention to the low-temperature case 
with T <C ep = h?kp/2m. Within this approximation, the discontinuous jump of the occupation 
number at the Fermi surface is slightly softened at low temperatures. Hence, the low-temperature 
expansion yields a power series in the dimensionless small parameter z = 7rT/2ei?, with expansion 
coefficients containing the derivatives of the zero-temperature result with respect to kp. Non-zero 
temperature introduces a new length scale, the thermal wavelength Xp = h-\/2'K/rriT. The small 
parameter of the low-temperature expansion can then be written as z = 2^“^/{kpXpy'■ In terms of 
At, we will be working in the regime of Xpkp ^ 1. 

The low-temperature expansion of the Green’s function, given by a one-loop Feynman diagram, 
produces artificial divergences where the spectral weight (96) is non-analytical. Fven though it is 
reasonable to assume that these singularities are not present when the full temperature dependence 
is taken into account, we continue here with the inspection of the Green’s function away from these 
singularities, only in the IR. The longitudinal and the transverse current-current Green’s functions 
take the following expanded form at low temperature. 


oo 

Go= ^ ao,j,k,m{ixy , 


(104) 


where the index o may either take the value ^ or t, indicating whether the Green’s function is 
longitudinal or transverse. The constants, read off from Eqs. (G31), are again listed in Table I. 
The different signs of a^ooo and a^ooi; shown in the Table, reflect the fact that thermal fluctuations 
tend to weaken the polarisation of the Fermi sphere and that the low-temperature expansion breaks 
down when the O term changes its sign. The inverse of the longitudinal Green’s function 

can be written in the form of Eq. (99), while the transverse Green’s function now takes the form 


with 


[GI]-' 


kph? 

mT'^ 


Ht 


(105) 


Ho 


OO 


j,k=0 


(106) 


where 

OO 

Kj,k = Y ( 107 ) 

m=0 


oo 

m=0 


The coefficients at low orders of the expansion are collected in Table IT The zero-temperature 
longitudinal Green’s function is invertible in the IR limit, enabling us to express its temperature 
dependence as a power series in z. Contrary to this situation, the presence of T 7^ 0 temperature 
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TABLE I. Coefficients of the retarded Green’s functions up to terms O O (y'^) and O (T^) in an ideal 


gas. 


jk 



^tj,k,0 


00 

-1 

1 

6 

0 

1 

8 

10 

TT 

2 

0 

TT 

4 

0 

20 

-1 

1 

2 

1 

13 

24 

01 

1 

1 

1 

1 

12 

24 

6 

288 

11 

0 

0 

TT 

16 

0 

21 

1 

5 

1 

1 

6 

12 

12 

24 

02 

1 

1 

1 

1 

240 

96 

60 

120 

12 

0 

0 

0 

0 

22 

3 

7 

7 

7 

80 

32 

240 

96 



TABLE 11. The same as 

Table I for the 

inverse retarded Green’s function. 


jk 

^i,j,k,0 

be,j,k,i 

bt,j,k,0 

^t,j,k,l 

^tj,k,2 

bt,j,k,3 

00 

-1 

1 

6 

8 

0 

0 

0 

10 

TT 

2 

TT 

6 

0 

— 1077 

0 

0 

20 

1-^ 

5 

6 8 

104 

3 

-64 

3277^ 

0 

01 

1 

12 

5 

72 

2 

9 

32 

3 

0 

0 


TT 

TT 


2877 

12877 



12 

12 

U 

9 

3 

u 

01 

1 TT^ 

13 77t^ 

124 

848 

512-1-40-n-’^ 

1 oc^2 


3 16 

18 96 

27 

9 

3 


no 

1 

1 

437 

224 

128 

0 


90 

45 

810 

135 

9 


12 

Ttt 

4377 

0 

26377 

12877 

25677 

480 

1440 


135 

45 

3 

22 

17 13-ir^ 

29 167^^ 

149 

1864 

7552 1 9747r^ 

1024 , 16647r’^ 

180 960 

72 5760 

810 

135 

45 ' 135 

3 ' 45 


is essential for the inversion of the transverse Green’s function, which can therefore be written as 
a power series in 1/z. 

We have now all of the necessary ingredients to return to the problem mentioned in Section 
VI A, namely, the understanding of the analytical properties of the spectral weight, needed to define 
simple collective modes. The point is that the spectral functions of an ideal gas are indeed analytic 
in the sufficiently small vicinity of x = y = 0. If the external source a^{p) is an analytic function 
and is negligible beyond the region of analyticity of the spectral weights, then the approximation 
of the Green’s functions, which is based on the analytical, small x and y behaviour, is justihed. 

The poles of the truncated Green’s functions clearly depend on the level of truncation. The 
truncation dependence is stronger in the transverse sector where the inversion of the Green’s 
function is possible only at finite temperature. The normal mode that arises from the O (x) 
truncation is 


x = -^z^ + 0{y^), 


(108) 
















35 


Im[x] 



Im[x] 



FIG. 10. The dispersion relation of the normal modes on the complex x = hmuj+{q)/kF\q\ plane, (a): The 
longitudinal mode for 0 < y < 0.4. (b): The transverse mode for 0 < ?/ < 0.2. 


which is to be contrasted with the O result, 

SiTT ± \/277r2 - 72 + 39^^ 


X = 


24 - 127r2 - 13.2^ 




(109) 


In other words, the transverse normal mode is fully damped at the O (x) truncation as can be 
seen from (108). On the other hand, the O (x^) truncation gives rise to a damped, but also 
oscillatory normal mode. As for the truncation of the series with respect to y, we note that the 
qualitative behaviour of the normal modes is already stable at the O (y^) truncation. The O (x^) 
dispersion relation qualitatively follows the numerical curve and is quantitatively wrong by a factor 
of approximately 2.5 for 0 < y < 1. The truncations at the level of O (x^) and O (y^), used for 
simplicity in all of the Tables and Figures, thus already provide a qualitatively correct description 
of the dispersion relations. Since |3f?[x]| < 1 at low temperature, the transverse sound wave is within 
the hydrodynamical regime. If we had used higher orders of x, the functions (C20) and (C23) could 
be approximated with even better (quantitative) precision and the plots of the truncated Laurent 
series would lie closer to those of the exact linearised equations of motion. 

The real and imaginary parts of the normal mode frequencies, u±{q) = ±ujr{q) + iu}i{q), are 
plotted in Figs. 10 for z = 0.87. The real parts were found to follow a quadratic dependence in the 
wave vector, uji^\q) ~ uJo£ — hq‘^/2m^ and ijj^\q) ~ cjot + respectively, with m|, > 0. 

An upper bound on the radius of convergence for the analytic form of Eqs. (99) and (105) is given 
by |a;+(q)|. The collective mode (109) is therefore a strongly damped, collective sound wave with 
both the speed of propagation and the inverse lifetime of the order O (T^). The damping of the 
normal modes makes condition (II), mentioned at the beginning of this Section, satisfied in the 
effective theory. 


D. Excitations in the presence of interactions 

In this section, we will include the photon-fermion one-loop electromagnetic interactions of the 
type depicted in Fig. 3 into the calculation to show that not all interactions are able to change 
the qualitative behaviour of the non-interacting gas. This can be achieved by modifying the linear 
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operator £ from Eq. (83) to contain the Coulomb interactions in its longitudinal component, while 
the contributions of the interactions to the transverse sector are negligible in the non-relativistic 
limit. We can further improve the precision of the one-loop calculation by making use of the full 
one-loop re-summed photon propagator, given by Eq. (76) and represented diagrammatically in 
Fig. 5. In this work, we will not include the vertex corrections nor the electron self-energy to the 
current-current two-point function, which would qualitatively change the behaviour of the electron 
gas and bring its dynamics closer to the usually discussed regime. 

The Schwinger-Dyson re-summation results in the replacement of by 

£, = - 


g2 - e'^Gl ■ 


The Laurent expansion now yields 


oo 


A = 2^ h,j,k,m,n{ixyy aokrY 

jkmn=Y) 


( 110 ) 


( 111 ) 


showing the emergence of a new length scale related to the Coulomb interactions. The coefficients 
bi,j,k,m,n are listed in Table III. We note that the partial Schwinger-Dyson re-summation of the 
photon propagator produces a non-perturbative result at finite temperature and density. Further¬ 
more, the longitudinal equation of motion changes at the C)(q*^)-order due to the expression (110). 
It changes the equation by taking —)• 2Ci. Intriguingly, the Coulomb interactions only intro¬ 
duce minor quantitative changes to the equations of motion and the sound wave remains the only 
collective mode at long wavelengths. The speed of sound is higher than in the ideal gas while the 
damping is only weakly influenced by the instantaneous Coulomb interactions, cf. Fig. 11. 

The inclusion of other one-loop corrections, the self-energy of electrons and the vertex cor¬ 
rections, would qualitatively change the hydrodynamical regime. They would generate a finite 
life-time and render the infrared limit, cj —)> 0, \q\ —)• 0, of the loop integral (C13) well defined. 
The loop integral would include an improved electron propagator whose complex self-energy would 
regulate the l/|q| singularities, thereby restricting the ideal gas hydrodynamical regime to 


T^hydr — {(l/j^) | 11^1 2//^l ^ 1 Ay > j (H^) 


instead of the domain presented in Eq. (95). Here, r^fp denotes the mean free path. These types 
of interactions open the way to the traditional phenomenological approach, which is based on the 
hydrodynamical regime with y < l/kprmfp and an arbitrary x. 


E. Equation of motion in momentum space 


To write down the equation of motion, we use the deviation of the density from the homogeneous 
value, n{q) = [j*^(g) — j'^(O)] /c, split the current into longitudinal and transverse parts, j = je+jt, 
where V • ji = 0, giving us 


VF 
Ti'^h ‘ 


j,k=0 


i,k 


lUJ 


VF\q\ 
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TABLE III. Coefficients of the longitudinal equation of motion of the Coulomb gas up to terms O {x^), 
O {y-^) and O (T^) . 
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FIG. II. The dispersion relation for the longitudinal mode in an interacting electron gas, plotted on the 
complex X = hmuj+{q)/kF\q\ plane. The range of y is the same as in Fig. 10, 0 < y < 0.4. 
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The temperature and the coupling constant dependence are included in the coefficients bijk and 
btjk- As a reminder, we are using the parametrisations = {(j)lc,a) and = {pc,j) of the 
external source and the vacuum expectation value of the Noether current. 

The continuity equation, dtn + V • ji = 0, can now be used to rewrite the equation of motion 
for the full current in the form of 


bt, 2,0 VF ujq 
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where the coefficients bj,k are given by 

k ^b2,o, , 

Oj,k = T - be,j,k - bt,j,k- 

O£,2,0 



(115) 

(116) 


The bj^k have been chosen so as to have a simple O (w^) piece. As before, the q ® q notation 
indicates the Kronecker product, giving us a matrix acting on j. 
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It is important to note that the linearised equation of motion breaks time reversal invariance. 
The pieces with even and odd powers of iuj arise from the and 51?S'2, respectively, in terms of the 
classification introduced in Section IVB, cf. Eq. (85). The hnite strength of in the linearised 
equation of motion generates true irreversibility because is odd under time reversal. Its presence 
leads to damping and relaxation to equilibrium, which are reflected in the non-conservation of the 
energy-momentum tensor, associated with the gas [9]. 

It is instructive to look for the stationary transverse flow with a; oc x = 0 and = {0, j{qy),0, 0), 
which is driven by = (0, a{qy), 0,0). It satisfies the equation of motion 


with the solution 


T2 

vph 


o-iq) — ( bt,ofl - \bt, oMjy + bt,o,2jj- ] j{q), 


kp 


‘k% 


(117) 


T‘^ a{q) 

VFh{q^ - q\){q^ - qV)' 


(118) 


where 


2 7 2 

— kp 


I 0,11 

26i,o,2 



(119) 


The argument of the square root is negative and therefore the stationary flow is a Gaussian wave 
packet in coordinate space with the width Ax = l/l^gil. Note that the role of the driving force 
of the flow, the pressure in hydrodynamics, is played by the external source in our scheme. 
The real and the imaginary parts of the wave vector are shown in units of fei? as functions of the 
temperature in Fig. 12. They characterise the characteristic length scale of the flow pattern and 
the dissipation of the flow, respectively. The lesson here is that the inverses of both length scales 
increase linearly with the temperature when the thermal energy is well below the Fermi energy. 


F. Equation of motion in space-time 


It is easy to find the equations of motion in space-time. All we need is to work out the operators 
that correspond to the multiplicative factors = |q|”, acting on the coordinate-dependent 

functions. Straightforward integration yields 



( 120 ) 


( 121 ) 


( 122 ) 


The short distance singularities of the Fourier transformation can be regulated by adding an in¬ 
finitesimal imaginary piece to the radius, which amounts to the use of the principal value prescrip¬ 
tion in carrying out the integrals. 
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FIG. 12. The location of the poles of the transverse current in Eq. (118) on the complex y plane, at a: = 0, 
plotted as the functions of the temperature T, in units of m = 1 and kp = 3.8 x 10“^. The colouring 
indicates the values of the temperature that grows from blue to red colour. At the origin of the y plane, 
T = 0 and we use the blue colour. The endpoint red colour is at T = 10“®. 


The equations of motion, written as 

= 

‘ft 

cf. Eqs. (91) and (92), have 
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(126) 


The linear response formulae for the components of the current, generated by an external 
electromagnetic field, which normally lead to the Kubo formulae, give in terms of a^. Our 
linearised equations of motion, i.e. Eq. (123), are similar to the Kubo-type equations, except 
that we now have expressed in terms of J^. The inversion of the kernel of the linear response 
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formulae provides an additional piece of information, beyond the Kubo formulae, namely the 
linearised equation of motion for J^, which is needed to make contact with the hydrodynamical 
evolution equations. As long as this equation can be considered local in time, it characterises the 
dynamics of the electron gas alone for t > 0, that is after the external source, o^, has been switched 
off, as explained in Section IV G. 

The kernels of the linearised equations of motion are given up to O (x^) in Eqs. (124) and 
(125). When higher orders are retained, these kernels contain non-local terms. In fact, the O (x^) 
contributions to the kernels are given in terms of the Fourier transforms of |q|"' , with n' > —n. This 
results in the appearance of the smearing kernels Q~^'■, which are of the order O (^\x — , 

for n' —3, and proportional to In \x — y\, for n' = —3, cf. Eqs. (120)-(122). The matrix elements 
of the operators Q~'^', with n' < —2, grow with spatial separation at a fixed frequency. However, 
such a non-locality is suppressed by the frequency in Vhydr, see Eq. (95), where the matrix elements 
remain bounded. Systematical and convergent improvement of the equation of motion, within the 
hydrodynamical regime, can thus be achieved by including higher powers of x. 

As an aside, we note that the long-range correlations, experienced by the composite sound 
modes, are reminiscent of the gluon dynamics in the confining vacuum of Yang-Mills theories 
where the single gluon modes develop strong correlations beyond a characteristic distance scale 
and glueballs are formed. Note that Q~^{x,y) grows linearly with \x — y\ so that the composite 
sound is localised with a linear potential at this order of the Laurent series. Higher orders generate 
correlations that increase even faster with spatial separation. 

Returning to the dynamics of the electron gas, the boost invariance requires special care because 
the truncation of the expansion of the effective action in the fields breaks the Galilean boost 
symmetry. A minimal extension that can be understood as a partial re-summation of the functional 
Taylor expansion to reinforce the Galilean boost-invariance, is provided by the replacement of the 
time derivative by a “covariant”, convective derivative dt ^ dt + v . It results in the replacement 

dt^ dt + (127) 

Po + n 

in our case with po = kp/Sir'^, and makes the equations of motion non-linear. The novel feature of 
these equations, compared to phenomenological hydrodynamics, is the presence of the factors 
and the two constants 6£,o,o and 6t,o,o- 

It is interesting to again, as in Section VIE, consider a stationary flow, this time in position 
space. The linearised equation of motion operator for this type of a flow takes the form 
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(128) 


The appearance of the corresponding operator of the Navier-Stokes equation. 


4s,.„ = -'f'5”A-(c+5) v■v^ 


(129) 


shows the presence of shear viscosity and bulk viscosity terms, which are smeared in space. This 
is of course not surprising from the point of view of an effective field theory. 
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Finally, we can look for the flow corresponding to a homogeneous external field, a' = (a, 0,0), 
(/) = 0, assuming the form of j = {j{y),0, 0) and n = 0. The solution of the equation of motion is 


a 

btfl,o 


f7=di 


(130) 


It contains two undetermined real parameters and is valid at any order of the truncation in 
y. The constants q± were defined in Eq. (119). The necessity to retain at least the O (y^) terms 
in the above equation of motion can be seen by noting that the O (y^) equation is solved by 


J = 


yt,o,o 


+ j+ cos 



+ j- sin 



(131) 


The current in Eq. (131) is purely oscillating. The level of the truncation thus clearly influences 
the position of the poles of the Green’s function and changes both qualitative and quantitative 
features of the solution. While quadratic truncation leads to a reversible oscillating current in 
(131), the quartic truncation introduces an imaginary part into the poles (cf. Fig. 12) leading to 
damping. 


G. Decoherence and irreversibility 

The build-up of decoherence in space-time is governed by the imaginary parts of the spectra of 
the quadratic effective action, isgiii{io)QC£{uj,q) and —f sgn(a;)9‘£t(u;, q). The imaginary parts of 
the spectra of Eqs. (C41) and (110), plotted in Figs. 13 and 14, vanish where the particle-hole 
spectral weights vanish, as demanded by Eq. (35). The longitudinal and the transverse decoherence 
of the density of the ideal gas at vanishing temperature, depicted in Figs. 13 (a) and 14 (a), show 
that the longitudinal density and the transverse current modes are classical at short and long 
distances, respectively. It is remarkable that the effect of temperature is qualitatively different in 
the longitudinal and the transverse sectors. In fact, the thermal fluctuations make the decoherence 
of the longitudinal modes stronger as one can see by comparing Figs. 13 (b) and 13 (c). The deep 
valley on the latter plot is due to the artificial singularity of the low-temperature expansion at 
|r| = 1 and should be ignored. The decoherence of the transverse modes becomes stronger with 
the temperature at higher frequencies, but it is strongly weakened in the infrared, according to 
Figs. 14 (b) and 14 (c). Note that the composite sound experiences strong decoherence both in 
the longitudinal and the transverse sectors. 

Decoherence and the breakdown of time reversal invariance share the same dynamical origin 
in a quadratic theory, as was pointed out in Section IV G. Hence, the effective dynamics of the 
current experiences decoherence and breaks time reversal invariance even in an ideal gas. The 
latter symmetry breaking manifests itself as the damping of the normal modes. When the system 
is placed into a finite space-time volume, then its spectrum becomes discrete and the imaginary 
part of the current-current Green’s function (Gl) becomes infinitesimal, 0{e). As a result, the 
far field component of the quadratic part of the effective action for the current vanishes according 
to Eq. (36) and the dynamics of the current becomes coherent and reversible. The emergence 
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FIG. 13. The longitudinal decoherence strengths, isgn{uj)QCi{x,y), as functions of x and y in the following 
cases: (a): T = 0, e = 0, (b): zoom into (a), (c): T > 0, e = 0, (d): T > 0, e > 0. Note that the singular 
lines in (c) and (d) are caused by the derivatives appearing in the low-temperature expansion in Eq. (C28). 


of irreversibility in the thermodynamical limit is another piece of circumstantial evidence of the 
spontaneously broken time reversal invariance. 


VII. DISCUSSION AND CONCLUSION 

In this work, we studied the equation of motion for the electric current in the low-energy, non- 
relativistic regime of an electron gas at non-zero density. Contrary to expectations, we found that 
the dynamical properties of the current are highly non-trivial even in the non-interacting, ideal 
gas limit. The simplest way to understand this behaviour was presented at the beginning of the 
paper where we argued that non-linear transformations of coordinates generate inertial forces, thus 
making the dynamics of non-linear coordinates (or fields) complicated and far from obvious. 

We showed that the electric current, which is bi-linear in the fundamental free electron field, 
obeys a highly non-linear equation of motion even in the absence of electromagnetic interactions. 
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FIG. 14. The plots of the transverse decoherence strength, z sgn(w)5£t(x, y), as functions of x and y in the 
following cases: (a): T = 0, (b): zoom into (a), (c): T > 0. As in Fig. 13 (c) and (d), the singular line at 
non-zero T is again caused by the derivatives appearing in the low-temperature expansion in Eq. (C28). 


The linearised equation of motion for the electric current was first derived in Fourier space. We 
found it to include dissipative terms and to display a non-analytic structure at vanishing frequency 
and momentum. The collective modes were sound waves, both in the longitudinal and the trans¬ 
verse sectors. It is important to stress that these modes are neither zero nor first sound modes of 
the usual Fermi gas. This composite sound mode, which was to our knowledge thus far unknown, is 
a result of the non-trivial redistribution of the energy-momentum, injected into the system through 
an external perturbation of the relevant composite operator, i.e. the Noether current. Our lin¬ 
earised equation of motion, an inhomogeneous linear equation for the current, is the inverse of 
the Kubo linear response formula for the current in the presence of an external vector potential. 
Our results are therefore equivalent to those derived from the linear response formalism. There is 
however an obvious difference between the two ways the dynamics is presented in the two schemes, 
namely, the equation of motion displays the collective modes in a more explicit manner than the 
linear response formulae. This fact may explain why the modes identified in the paper had not 
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been found previously within the usual linear response formalism. Decoherence, which goes beyond 
the linear response formalism was also analysed in this work. We found it to be to be very strong 
for the composite sound modes, thus making them important also in the classical limit. 

The equation of motion for the longitudinal component of the current was found to be local in 
time for an ideal gas at vanishing temperature, so long as the motion was slower than the Fermi 
velocity. We could recover a similar local equation for the transverse current, but only at non-zero 
temperature. However, both the longitudinal and the transverse equations of motion remained 
non-local in space. 

After obtaining the linearised equation of motion for the current, we studied its non-linear exten¬ 
sion by replacing the time derivatives with convective derivatives. This minimal extension enabled 
us to restore Galilean invariance. Although the equation of motion for a stationary flow showed re¬ 
semblance with the Navier-Stokes equation, important differences were present in comparison with 
phenomenological hydrodynamics. In particular, our system of “mechanical” equations was closed 
without evoking any thermodynamical considerations. Furthermore, spatial non-locality arose due 
to the factor of l/|q|. This non-locality followed from the singularity of the current-current two- 
point Lindhard function at cu = g = 0, generated by the gapless particle-hole excitations. 

The effects of the microscopic QED interactions on the dynamics of the current can be most 
easily assessed in the framework of the skeleton expansion. Intriguingly, what we found was that the 
re-summation of the Coulomb interaction into the photon self-energy only changed the linearised 
equation of motion quantitatively, while preserving all of the qualitative features of the behaviour 
of a non-interacting electron gas. We therefore determined that qualitatively different behaviour 
would have been generated by vertex corrections and the electron self-energy insertions. Thus, 
such corrections would be required to bridge the gap between the electron gas regime studied in 
this work and the interacting regime that is usually discussed in physics. 

At the end of the paper, we turned our attention to the discussion of decoherence and ir¬ 
reversibility. We argued that it is the non-vanishing far field component of the current-current 
Green’s function that connects irreversibility with decoherence. This relation could be interpreted 
as a generalisation of the fluctuation-dissipation theorem. As a result, the longitudinal and the 
transverse modes displayed strong irreversibility and decoherence at high and low frequencies, 
respectively, even in the absence of interactions. 

We believe that beyond our presentation of the intricacies and complexities that govern the 
dynamics of non-interacting, ideal gases, the results presented in this work have methodological 
value and may lead to the development of a systematic way to derive hydrodynamical equations 
for more involved models. Furthermore, our calculation gave us the exact linearised equation of 
motion with irreversibility and decoherence; a result that would not be invalidated by calculations 
at higher orders. 

This work opens several questions: What is the signature of the free, non-interacting nature 
of a dynamical system when expressed in terms of non-linear functions of the coordinates? The 
structure of the connected Green’s functions completely disguises the underlying harmonic system. 
Is there perhaps a simpler way to recognise the existence of the uncorrelated degrees of freedom than 
by studying the properties of its Noether current? How precisely do the interactions, in particular 
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the vertex corrections and the electron self-energy change the picture? Does our approach allow 
us to recover the full analyticity of the equation of motion for a Fermi liquid in Fourier space, 
assumed in phenomenological approaches? 

Beyond these issues, what we find to be most important at present is to assess the relevance of the 
phenomena discussed in this work for the actual experimental observations. The central question 
is whether the correlations, established by the ideal gas equation of motion, remain experimentally 
accessible in an interacting system. To address this issue, one should distinguish between three 
independent length scales. The first is the average inter-particle separation, Xjkp. The interactions 
introduce the second length scale, the mean free path, r^fp- The experimental observations are 
then carried out at the third length scale, Vq. The three physical scales normally obey the relation 
1/kp < Vmfp < Vo- The elementary particle regime is restricted to the length scales i < 1/kp. In 
this regime, the energy-momentum is mainly carried by the particles. This is because the holes, 
which have \q\ < kp, carry a relatively small fraction of the total energy-momentum. Therefore, 
no collective modes can exist for £ < 1/kp. The ideal gas collective modes appear at £ > 1/kp, 
but they remain restricted to the composite length scale regime, 1/kp < £ < r^fp, because the 
particle collisions generate a finite life-time for the dressed quasi-particle modes for i > r^fp- The 
usual experimental observations at Tq allow one to reconstruct the physics of the quasi-particle 
regime, r^fp < £ < Tq, which is dominated by the “true” quasi-particles of the Fermi liquid. Their 
hnite life-time screenes the l/\q\ singularity of the ideal gas equation of motion, thus justifying 
the usual phenomenological treatment. The composite sound waves are therefore the collective 
modes that play a central role in the composite particle regime by forming the internal structure 
of the quasi-particles. However, they are suppressed on scales beyond the mean free path, i.e. in 
the quasi-particle regime.^ To uncover the physics at the composite particle scale, we would have 
to resolve the structure of the quasi-particles and the effective couplings. The main question is 
therefore whether there is an experimental method to zoom into the short distance features of the 
dynamics in order to identify the collective modes, down to the microscopic edge of the collective 
phenomena, at 1/kp. We leave the resolution to this question open in this work. 

We further believe that it would be interesting to extend the calculation to the relativistic 
domain where dynamical degrees of freedom of the electromagnetic field appear and one could 
systematically derive equations, which are analogous to magnetohydrodynamics. We expect that 
the radiation reaction would enhance the dissipative forces and generate new collective modes. 

Finally, we may also wonder if turbulence can occur in our system when the equations of 
motion include the minimal non-linear extension, cf. Eq. (127). As far as dimensional analysis is 
concerned, the characteristic scales of the non-interacting electron gas at zero temperature are the 
mass m and the average particle separation, 1/kp. We can thus construct a Reynolds number, 

ukp£ uh£ 

l\p 5 

vf Tn 

® We find it interesting that there exists a similarity between the composite sound and the dynamics of gluons in 
QCD. This is because both the single gluon excitations inside the glueballs, i.e. the quasi-particles of QCD, and 
the composite sound modes are localised into a finite region in space by strong correlations that increase with the 
distance. Outside the quasi-particles, the modes are screened. In QCD, this is due to confinement. 
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where I and u denote the typical length and velocity scales of the flow. One finds that Rp ^ ui in 
CGS units. By assuming the flow velocity to be of the order of the magnitude of the metallic Fermi 
velocity, u = vp ^ 10®cm/s, we find that Rp ~ 10®£. In viscous fluids, turbulence sets in at around 
Rp = 4000. If such a naive dimensional analysis can be carried over to our case then one expects 
that a homogeneous flow of the size i ~ 400nm indeed reaches Rp ^ 4000 and thereby experiences 
turbulence. The presence of a heat bath introduces further length scales, the thermal wavelength 
\p and the Coulomb interaction brings in the Bohr radius, oq, rendering the dimensional analysis 
more involved. In fact, these gases have two further Reynolds numbers, Rp = Rp/{XTkp)‘^ and 
Rc = Rp/(aokp)'^, where Rp characterises the ideal gas at finite temperature and Rc the Coulomb 
gas. However, Rp <C Rp in the low-temperature limit, Xp 3> l/kp, and Rc < Rp for a weakly 
coupled Coulomb gas, where oq 1/kp. Hence, Rp remains the relevant Reynolds number for 
determining the onset of turbulence. 

We end this paper by noting that it would be fascinating to connect our results with the low- 
energy analytic structure of a field theory dual to a higher-spin Vasiliev theory. The reason for 
this possible analogy is the fact that the presence of higher-spin conserved currents constrains the 
field theory correlators to be those of a free theory of bosons or fermions [43] . We hope to explore 
some of the research directions and open questions raised in this work in the future, particularly 
in models with direct relevance for experimental observation. 
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Appendix A: Linear response and hydrodynamics 

The linear response formalism [21] is a widely used technique for calculating the expectation 
value of an observable, Op, in the Schrodinger representation, provided that the Hamiltonian is of 
the form H = Hq -|- Hi{t). The expectation value, 

{{Om = Trp(t)0(t), (AI) 

can be written in the interaction picture representation, using Hq and Hi{t) as the non-perturbed 
Hamiltonian and its perturbation, respectively. In this representation, an operator A is then given 
by the expression 


Ai{t) = efi^°Ase 


(A2) 
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and 


Pi{t) = U{t, ti)poU\t, ti) (A3) 

is the density matrix, where the time evolution operator, satisfies the equation of motion 


ihdtU{t,ti) = 

The equation of motion for the density matrix, 

ihdtpi{t) = [Hu{t), pi{t)], 

can be solved iteratively, 

Pi{t) = Po-T [ dt'[Hu{t'),po] + O {Hfi) . 
J ti 

The first, linear approximation in Hi gives 

= mm - T,pmt) 


(A4) 


(AS) 


(A6) 
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dt'Tvpii[Hii{t'):Oi{t)]. 


(A7) 


One usually assumes the separability of the time dependence in the perturbation, Hi{t) = 
and writes 


smt))) = - rdt'G^mmt'), 

Jto 

where the retarded Green’s function is defined by 


(AS) 


(A9) 


To relate this result to the CTP formalism, let us suppose that we want to find the expectation 
value of a local composite operator, F{(p{x)), in the scalar field theory, defined by the action 


Sh[(p] = S[(J)] + / d‘^xa{x)G{(J){x)). 


(AlO) 


The perturbation series of the expectation value of our composite operator in h can be calculated 
by means of the generator functional 


as 
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(A12) 



48 


for either a = + or a = —. The linear response formnla is the O (h) result, 

,52 52 


6 {{F{cPix))) = j d^y 
which can be written as 


JW[a,j] 


5 f{x) 5 a+{y) 5 j^{x) 5 a {y) ^ 

5 {{F{(t){x))) = f d‘^yD^{x,y)a{y), 


|a+=—a“ =0.^=0’ 


where 


nr/ / ^‘^W[d,j] 

5j^(a:^)5a'"'(y)|a+=-a-=a,i=o' 


(A13) 


(A14) 


(A15) 


Hydrodynamics addresses the inverse problem. There, we are interested in the equations (of 
motion), satisfied by the expectation values, in which the external sources appear linearly. In case 
of the linear response theory, it is easy to find the equation in question. 


a{x) = - I d^y[D^{x,y)]-U{{F{cjriym. 


(A16) 


The only subtlety is the potential necessity to exclude the null-space from the domain of the inverse 
Green’s function. 

The generalisation of such an inverse linear response formula beyond O (a) is provided by the 
functional Legendre transform of IT[a, j], the effective action. 


r[F] = W[d,J]-aF, 


(A17) 


where 

6 F 

The inverse Legendre transform is then also given by Eq. (A17), with 

mF] 

a =-^—. 

dF 

This equation plays the role of the equation of motion and produces a non-linear extension of the 
hydrodynamical equations. The inverse Legendre transform generates the non-linearity, which is 
necessary to close the equations, without an introduction of auxiliary variables, such as thermody¬ 
namical functions. 
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(A18) 


Appendix B: Free propagators 


The detailed calculation of the free propagator is easiest to carry out in the conventional operator 
formalism, see e.g. [28, 29]. We summarise the results below. The boson propagator, given by Eq. 
(26), yields 


D{k) 


( 


—2'Kid{k^ — 




+i€ 

2^2 


-27rid(k^ - ^)e(-k^) 


~w 


)0(A:°) 


— i2Tr5 ( k 


fc2_: 




J 


nsik) 



(Bl) 
















49 


in Fourier space where the Bose-Einstein distribution function is 

0(-fc°) 0(fc°) 

— 1 — l" 


The inversion, cf. Eqs. (34) and (35), gives the following expressions for A = Z) ^ 


A”'(Z) = A:^ — m^, A*(/c) = e, A'^(fc) = z sgn(A:°)e. 


(B2) 


(B3) 


The free fermionic propagator, defined by the generating functional 


et^b'4l 




(B4) 


can be written as 


G^^{x,y) 


^{0\T[^|J»{x)^P<^{y)]\0) -{0\^pf^iy)i^<^{x)\0) \ 

[ {oirixWiym {o\T[{^^Hy))HHxh^rm*) ’ 


(B5) 


and the detailed expression written in terms of the scalar propagator in Eourier space is 

G(fc) = 

In case of finite temperature and density, one uses the occupation number density, 

0(fc°) 0(-fc°) 

' g/3(efc-M) _)_ 1 e/3(€fe+zi) _)_ 1 ’ 

and the full propagator becomes 


(B6) 


(B7) 


Gk 


(^ + m) 


Dk + 27ri6{k'^ 


m‘^)nF{k) 



(B8) 


Appendix C: Current-current two-point function at finite density in the non-relativistic limit 


The result of the calculation of the current-current Green’s function, given by Eq. (78), 

acW) = -' / (0i‘' . (ci) 

at finite density and vanishing temperature in the non-relativistic limit, c —)■ oo, is briefly sum¬ 
marised in this Appendix. The real and imaginary parts if the Green’s function are always defined 
in position space; hence, 

^xGiq) = ^ [Giq) + G*i-q )], z9,G(g) = ^ [G(g) - G*(-g)]. 


(G2) 
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1. Lorentz structure 


The two-point function is symmetric, = G(o-'i/)(o-/^)(“P)) and transverse, p^G^^ip) = 

0. Furthermore, it is covariant and depends on two four-vectors, p^ and /3^, defined as (3^ = (1, 0) 
in the inertial frame where the electron gas is at rest. This gives two independent kinematical, 
scalar combinations, = — [q — u{uq)]‘^ and ^ = uq/\q\ = Lo/c\q\, where the notation = {u}/c,q) 
is used. Such a tensor can be parametrised by two Lorentz scalars as 

+ DtPr, (C3) 

where Pt and Pi are projectors onto the three-dimensional transverse and longitudinal subspaces, 

^0 o\ 


pr = - 


\^TJ 
1 ( 1 


pfip _ 

^ i-eweL)' 


(C4) 


respectively, with n = fc/|fc|,T = n(8)n and T = 1 — L. The inverse, dehned by G^pG 
is given by 




For future reference, the retarded and advanced photon propagators are 

1 1 


D^q) = -Dl,{q) = - 


{q^ + ie^ — q^ ’ 


DIM) = -DIM) = - 


{q^ — ze)2 — ' 


(C5) 


(C6) 


2. Vacuum contribution 


The two-point function is the sum of the vacuum and the finite density contributions, 

^ ^vac ' ^gas") 


(C7) 


and both the vacuum and the finite density contributions are of the form (28). The vacuum 
contributions, Gvac = GyacT^'^ ^ are easy to find. The diagonal CTP block, G/M — 
gives the standard result. 


G++(<?)-^g'|^ + 2(l + 


2m?c^ 


q^ 


1 4m^c^ 


— 1 arccot 


' 4m^c‘^ 


- 1-1 






g 

IStt 




+ 0 


.2 \ 2 


m?c^ 


(C8) 


The off-diagonal CTP block, calculated by using the free propagator (64), is G^^^^ = Gfy^^y = GM^ 
with 

G^~{g) = ^ j ^^^27r(5((p-hg)^-m2)0(-p°-g°)27r(5(/-m2)0(p°)trV(p-Fg,g), (C9) 

















51 


where the trace is taken over the Lorentz indices of the trace formula, 

q) = tT^'^ + ?7ic)7^(g/37^ + me) 

= - pq)g^'' + ^p^q’^ + 


(CIO) 


Simple steps lead to the expression 


= O I I ^" Q(-g - m) 

SttcIqI 7o 


2 ;^, 0 dpp 

t-si _ m) / 

ijj< 


— Q{2p\q\-\q^+2ujpq^\), (Cll) 


which can be neglected in the non-relativistic limit because the Heaviside function vanishes for 
non-relativistic frequencies with <C m. 


3. Fermi sphere contribution 

To find the non-zero density contributions to Gi and Gt in the electron gas we need to compute 
the CTP blocks corresponding to the time components, T = Ggasy the spatial trace, S = Ggas- 
For this purpose, we need to use the trace factors, cf. Eq. (CIO), 


t = N^°{p + q,p)\p 2 =^ 2^2 =t' - 2{q^ + 2pq), 
s = W\p + q,p)^p 2 =^ 2^2 = s' + 2{q^ + 2pq), 


where t' = 8{{p^)^ + p^q^) + 2q^ and s' = 8{p^q^ + PG — in the loop integrals 
, -l-l- 


(IS)^^-/«) 


27r5{{q + p) — c^)nq+p2'n5{p — m^c^)np 


(Uq) 

V‘S(9) 


. 2'k5{p^ — m?c^)np ,2Tr6((p + q)^ — mremup+g 

— % - — — — % - — — — 

{p + ~ rn?c^ + p^ — + ie 

=—I I [ ii) 27rd{{q + p)‘^ — m‘^c‘^)27T6{p‘^ — 


X [ 0 (-p° - q^)np + 0 (p°)np+g - Ug+pUp] . 

At this point, it is advantageous to introduce the integrals 

d'^p 


h[q-,f] = 
h[q]f] = 

and write Eqs. (C13) as 

/T\+- 


(27r)4 

d'^p 


a) 


f{p,q)2'K5{p‘^ 

f{p, q)2TT5{p^ - m^c^)2TT5{q^ + 2pq), 


Ur, 


q^ + 2pq + ie 


+ -/(2) 

2 


q'i [s' ) '^p^q+p 


+ {q^ -q), 


(s) (q) = q; (s ') (©(-/ - q^hp + Q{p°)np+q - uq+pUp) 


(C12) 


(C13) 


(C14) 


(C15) 
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These integrals are then evaluated in the non-relativistic limit, \p\ <C me, by assuming that the 
integrands are spherically symmetric and non-vanishing for > 0. We find 


/(i) 


9 ; 


g{p°,p) 


1 


+ 2 pq + ie 


/( 2 ) 




q;h{ {p^,p), {-,q 


dvr^mc ^ 

IGvr^lglmc Jq 
1 f 


dpp‘^f{mc^,p), 


dppg{m(?,p) log 


k + p + ie 
k — p + ie' 


IGyr^lglmc J v' ' c 

where the following dimensionless parameter has been introduced, 

+ 2mu} 2m(jj — 


d^ph {{mc^,p), ^{rkp -Pz), (CIG) 


r = 


2\q\kF 


2\q\kF 


(C17) 


We now need to find the Fourier transform of the real part of the CTP diagonal block at the 
leading order in 1/c, 



mi 


a) 


Ur, 


+ 2 pq + ie 


-|- (cj —)■ —oj) 


1 


IGvr^lqlmc , 


dpprip 


( Srn^c^ 

1 8(maj-|-p2)+2(jr2 


log 


k + p 


k — p 


+ 8|g|p( /) 


+ (w —)> —uj). (C18) 


The momentum integrals can easily be carried out at vanishing temperature, where Up = 
Q{p^)Q{kF — IpI), leading to the result 


mT++{q) 

3?x5++((?) 


k%mc ^ , 


k\ 


27r^mc|g| 


kpLzi'f') + 


moj + — ) Li{r) 


+ 


Gvr^m 


+ (w —>■ ~w). 


where we have defined 


Li{r) 

Lsir) 


dkk log 


r + k 


r — k 


= r+-{l- r^)log 


r -t- 1 


r — 1 


dkk ^log 


r + k 

r 1 ,, 4,, 

r -|- 1 

r — k 

= g + y + |(1 )log 

r — 1 


(C19) 


(C20) 


The next step is to compute the off-diagonal H— CTP contributions to T and S. At the leading 
order in the 1/c expansion, the off-diagonal CTP block is 

, +- 


(T{q) 

\S{q) 


= -il^‘^^ 

ime 
2Tr‘^\q\ 


[ 0 (p°) - Up] Uq+p 

d^p 


^(ma;+p2 + ^)y 


^ Q{kF -\p + q|) 0 (|p| - kF) 6 {pz - rkF). (C21) 


We choose the ^-direction to be parallel to q. This integral is over the region of a plane of the 
height pz = kFV, which is shown in Figs. 15 and 16. The integration thus needs to be taken 
over the points of the plane which are outside of the Fermi sphere, centred at the origin, and 
inside of another Fermi sphere, which is centred at —q. It is advantageous to parametrise these 
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FIG. 15. Domain of integration in Eq. (C21) for non-overlapping spheres, case (a), indicated by the solid 
line. 



FIG. 16. The same as Fig. 15 except for overlapping spheres, cases (b) and (c). 


integrals by the dimensionless variables x = and ?/ = |^- There are three different functional 
forms for these integrals, corresponding to the three cases shown in Figs. 15 and 16: (a): y > 2, 
—y — 1 < r+ < —y + 1, (b): y < 2, —1 — y < r_ < —1, and (c): y < 2, 1 < r+ < —where 
±x — I in the non-relativistic limit. Straightforward integration yields 


_ q^±2mu} 
^ 2\q\kF 


where 




, . imck% ,, , 

■( 9 ) = -y^T/i(x,y), 


S+-{q) = - 


2vr|q| 


27rmc|q| 


kpM 3 {x, y) + Lom + kpr^ + — Mi(x, y) 


rp2 

Mi(x,y) = 2 dpp= { 

Jpi 


M3{x,y) = 2 


rP2 I 

/ dpp^ = ^ 

I pi ^ 


1 — (a) 

1 — ( 6 ) , 

-2xy (c) 

{l-rlf 
{l-rlf 


(a) 

(5) . 


xy{y‘^ + 4x^ — 4) (c) 


(C22) 


(C23) 


The integrals over the interval are given in terms of p 2 = yjl — r‘^ and the relative locations of 

the two Fermi spheres are (a): pi = 0, (b): pi = 0, (c): pi = y^l — 
off-diagonal CTP block is therefore 


The real part of the 




fnq) 

\S{q) 


n-q) 

S{-q) 


(C24) 



















54 


4. Low-temperature expansion 


In this part of the Appendix, we calculate the loop integrals (C13) at a small, non-zero tem¬ 
perature, T <C Ci? = ^, to find the temperature-dependent corrections to T~^~^, and 

. The temperature dependence of the vacuum contribution is strongly suppressed and will be 
ignored. We start with the ++ CTP block, for which one needs to evaluate integrals of the type 

dxf{x) 


r 

I{T) = / 
Jo 


(C25) 

e T -^1 

The low-temperature expansion involves the primitive function F{x), where F'{x) = f{x), and is 
introduced because of the partial integration, 

F(0) 


I{T) = I{T) - 


-lEl 

e T + l 


/(r)-F(o), 


where the exponentially small quantities are neglected and 

poo 

i{T) = - / dxF{x)d^- 


e T 


roo 

- -= / dyF{^JTy + ep) 


{ey + 1)2' 


(C26) 


(C27) 


The next step is to expand F{y/Ty + ep) around y = 0, extend the integration over —oo < y < oo 
and keep the even part of the integrand, while the odd part vanishes. By neglecting the O (T^) 
terms, the leading-order result becomes 




2F{P^) + -[^— - 


To compute the remaining integral, we write 


y^ey 


v+ip’ 


(C28) 


(C29) 


and solve for n = 0 and n = 2, finding the values Kq = ^ and K 2 = By using the above 
expressions, we can find the leading-order result for the integral (C25), which is 


I{T) = 7(0) + [kldl^m) - kpdk.m] , 


T2 


(C30) 


showing that the low-temperature expansion is an expansion of the integrand around the Fermi 
surface. The final expression for the ++ CTP block, which includes an 0{T‘^) correction of the 
zero-temperature result presented in Eqs. (C19), is 


^ 4mc 
^ 27r2|q| 


3^.5++(g) = 


kl 


Li{r) + 


kpL^ir) + i mu + ^ ] Li{r) 


+ (w — > — w). 


27r2mc|q| 

(^F[8^3(0 - 5 r 4 (r) + r‘^L'^{r)] + (^mu + ^ ) [r^L'l{r) - rL[{r)] 

/ ^2 rp2 \ 


(C31) 
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The +— CTP block contains an integral of the type 


/(T) = / dxfix) 


e T 


ex +ieT +1 


(C32) 


which at low temperature receives non-zero contributions from the interval ei < x^ < € 2 , where 


ei = e_p — 


kpx^ 

2m 


€2 = cf — 


{qlkpr k‘pr‘^ 


2m 


m 


2m 


After partial integration, we find 


I{T) 


HT), ei < 0, 

ei>0, 


where 


/(T)= r dyF N{y), 


and 


N{y) = 


+ 


(C33) 


(C34) 


(C35) 


(C36) 


{l + e''i-y)'^{ey-'^2 + l) {ey-’'^+l){l + e'^2-y)2 (1 + 1) ’ 

with Uj = Y and v = The function N{y) in the integrand has a peak around y = uj] hence, 


N{y) 




(l+e'"!-!') 


(C37) 


(l+e'"2-!^)2 ) y ^ ^ 2 - 
By expanding F{y^Ty + ej), the integral becomes 


/ dyF{y/Ty+Vj)N{y+Uj) « / dy 
l-Ay Jo 


1 


^ ( F'\v^) - 


F' 


y2T2 


(1 eyy 

(C38) 


We take Ay to be large, but smaller than 1/T, and find 


/(T) = Fi^2) - F{^,) 
24 


2 

+T^ — 




1 




(C39) 


n +1 


Now, the expressions for Mn in Eqs. (C23) contain the functions Fn{x) = 27r(mx^)^“, with n = 1 
and n = 3, when written in the form of Eq. (C32), i.e. at finite temperature. By using Eq. (C39) 
and the fact that only n = 1 and n = 3 are relevant, we recover 


4(r) = FniVe2) - i4(VeI) + F 


TT 


1 


62 


- F"(V 62 ) - 




) ei 


- - I Kisixo - 
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= 4(0)+T' 
= 4(0). 


,TT'^n-l-l fn + 1 

y 2 


— 1 m 2 


n + 1 


-2 


--2 
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^4I 

(C40) 
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What this result implies for the low-temperature expansion of Eqs. (C22) is that the finite temper¬ 
ature dependence of the spectral functions and only appears at the O (T^) order. For 
the purposes of this calculation, we neglect the O (F^) contributions, which makes the temperature 
dependence of the -|— CTP blocks drop out altogether. 

Finally, the retarded Green’s function, Gr = Gn + Gf = G++ — G_|_, has longitudinal and 

transverse components, as in the parametrisation of Fq. (C3). In the non-relativistic limit, the 
two sectors are dehned by the expressions 

Gi = ^ 3?, [r++ - r+-], 

= 13^. [(r++ - r^-)e - + 5+“]. (C41) 

The variable ^ was defined above as ^ = w/c|q|. 
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